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PREFACE 

This book is based on a course given in recent years to 
students in Chemical Knginoering at the Massachusetts 
Institute of Technology. It is similar in grade of difficulty 
to the "Differential Equations 1 ' of Dr. II. I*. Phillips, from 
which it differs in containing more Chemistry and less 
Mechanics. 4 

It would be impossible 1 to mention all the sources from 
which the problems have been taken. A large number of 
them are due to Professor W. K. Lewis, who was the first 
to rocogni/c the need of such a course* and who has given 
the authors many valuable suggestions. Some of the simpler 
problems on rate of reaction have been adapted from "Chcm- 
icjil Statics and Dynamics" by J. W. Mellor, a work to bo 
commended as collateral reading. The 1 mrtlwd erf treatment, 
however, has for the iimst. part evolves! in the actual course 
of class-room teaching, nnel is belie 4 veel to be* largely new. 

The purpew of the be>e)k is not to teach Chemistry, 
but to teach Mathematics in a form readily assimilated 
by Chemists and Chemical Engineers. Anel just sis Chem- 
istry itself ciiiinot be we*ll unelerstewel except* as a branch of 
Science in a wider aspect, so the principles of Mathematical 
Chemistry are? often most clearly brought out by the? study 
of problems much broader in their application than any 
specific technical requirement. The practical necels erf the 
Chemical Engineer are, however, constantly ke*pt in sight. 
Much emphasis is placed on graphic mesthoels anel on numer- 
ical work. After careful consideration the answers to the 
problems have with few exce*ptie>ns been omitted. The 
authors believe that a chemist or an engineer, or a young 
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man about to become an engineer, should possess sufficient 
stability of mind to check his own calculations. 

A few of the more difficult problems lead to differential 
equations of a type whose solutions have never been worked 
out. These may be of interest to the advanced reader in 
suggesting topics for research. 

The authors will be grateful to anyone pointing out mis- 
prints or other errors. 

FRANK L. HITCHCOCK 
CLARK S. ROBINSON 
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DIFFERENTIAL EQUATIONS 
IN APPLIED CHEMISTRY 



CHAPTER I 
INTRODUCTION 

1. Why Chemists require Mathematics. Mathematics 
may be defined as " that science which treats of numbers 
and magnitude, or of whatever is capable of being numbered 
or measured."* Its business is to examine the facts within 
its scope, that is, all processes or quantities which can 
be numbered or measured, and so to classify and analyse 
them that they may be used with advantage and with de- 
creased labor. Such a definition is very broad. If we 
accept it as true it follows that we shall include physics, 
chemistry, and engineering within the realm of mathematics. 
s These branches are indeed often called the mathematical 
sciences in contradistinction both to the fine arts and to 
those sciences which arc purely descriptive or historical. 
It is a very significant fact, however, that at the present 
day all phases of human thinking tend to become more 
analytical and more exact, for such is the very spirit of a 
scientific age. If there are exceptions, they occur precisely 
where we should expect them, namely in some developments 
of the fine arts; but such subjects as economics, biology, 
geology, and even anthropology, where formerly a rigorous 
quantitative treatment was almost unknown, are being 
brought more and more within the viewpoint, and to some 

* Worcester's Cowpre/iercswe Dictionary. 

1 



2 DIFFERENTIAL EQUATIONS IN APPLIED CHEMISTRY 

degree within the actual working methods of Hie mathe- 
matician. 

It may then be taken as the function of mathematics to 
assist in the classification and analysis of the facts regarding 
physics, chemistry, and the like; and its success in growing 
from an abstract science into a broader field of usefulness is 
measured by what can be accomplished by its aid. 

2. Calculus as a Tool. Mathematics in the rather vague 
language of common use is often divided into " lower " 
and " higher." There is no sharp line of division. Cal- 
culus is more fundamental than Analytical Geometry or 
Trigonometry. It is not only a method of exact thinking 
but includes labor-saving devices and short cuts for obtain- 
ing results. It may possibly be true that, as has sometimes 
been said, there is no result obtainable with calculus which 
could not (conceivably) be obtained without it. But the 
labor involved would frequently be prohibitive, and the 
methods would be so complex that no human mind could 
find them out. 

The need for labor-saving methods of thought in chemistry 
and chemical engineering is greater now than ever before. 
The calculus is a labor-saving tool of the greatest importance. 
The chemist or engineer who fails to make use of it in the 
solution of problems is handicapping himself. In these 
days of strong competition and of the elevation of the fittest 
he who masters its use will find himself amply repaid. 

The first object of this book is to help students in chem- 
istry and chemical engineering to think more readily in 
terms of calculus. Most engineering schools give excellent 
courses in differential and integral calculus, usually during 
the first two years of the four years' course. This is neces- 
sary in order to enable the student to handle readily the 
fundamental subjects of physics, physical chemistry, applied 
mechanics, and so forth. However, when the more special- 
ized work of the fourth year is reached he is better ready to 
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appreciate the bearing of differential equations upon the 
material in hand. It is not our purpose to review in any 
systematic way his previous mathematical studies, although 
references to them will be frequent and occasionally detailed. 
Still less is it intended to teach him in this little book the 
technical problems of his profession with any complete- 
ness, even though the examples given him to work are 
selected with direct reference to actual chemical practice. 
Some of them are based on investigations carried on at 
the Massachusetts Institute of Technology in recent months. 
The purpose is rather to lead him into fresh fields of thought, 
where he may find a broader point of view, and whence he 
may at the same time regain a stronger hold of the funda- 
mentals and learn to apply them to his professional work, 
discovering the value to himself of this wonderful tool. 
Perfect mastery in technical problems comes only with 
long practice. This book may serve to acquaint him at 
least with the existence of the tool. 

3. The Scope of this Book. The point of view of the 
authors may be indicated by an allegory. Suppose a geolo- 
gist and an artist to be taking a stroll together through a 
varied landscape. Each sees the same hills and valleys, but 
to one observer his habits of thought tend to present these 
scenes to his mind as so many scientific questions, while 
to the other they are more likely to be material for his brush 
and easel. The purpose of this book will be attained in so 
far as the student acquires the habit of seeing the phenomena 
of life and of nature, including his professional work, as so 
many examples which might be set up in the language of 
differential equations. 

Accordingly the reader will find in the following pages 
various problems which may seem at first remote from 
chemistry, questions on the growth of populations, on 
colonies of bacteria, on the age of the earth, and on the rate 
of drying of a wet sheet hung in the wind. As his familiarity 
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with the new mode of thought increases, however, he will 
come to see all these as cases of " mass action", words which 
have a familiar sound to every chemist. He will be able 
to say more exactly just what he means by the mass law, 
and he will note that the logic by which, from a study of 
radioactive changes in rocks, their age is estimated, is very 
similar to that by which he can learn to stop a chemical 
process at precisely the moment to obtain the best yield. Or 
from the simple matter of the wet sheet, he will work toward 
the complex types of drying apparatus used by the chemical 
engineer. And he will note that laws of growth and decay, 
whether they govern the rise and fall of a race of men, or of 
an empire, or of a culture of bacteria, are plotted by curves 
of a form very similar to those governing standard types of 
chemical processes. The fact that cultures of bacteria are 
themselves of great use in some chemical operations is 
interesting and important, yet of minor significance in 
learning the method of thought. The increase in popula- 
tion in early New England times is a still better, because 
very simple, illustration of the mass law, in spite of the 
fact that the Puritans possessed few chemical industries. 

The two-fold object which most of the problems in this 
book possess, of providing concrete illustrations of funda- 
mental principles while at the same time working toward 
the actual necessities of chemical engineering, is to be spe- 
cially noted in connection with the series of problems on 
the flow of fluids. This series runs through the whole book. 
It begins in Chapter I with the problem of washing brine 
out of a tank. There does not seem at first much to be 
gained by wasting the good salt in this way. But in the 
first place these problems on flow offer the best mechanical 
model of " consecutive processes," witness the use made 
of them in that fascinating book " The Interpretation of 
Radium " by F. W. Soddy. And in the second place they 
introduce the student directly to chemical processes of 
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the continuous current type which are commercially im- 
portant. 

To sum up what has been said thus far, the first object 
of the student should be to learn the language of calculus, 
so that he may express all these problems in the form of differ- 
ential equations with the same ease with which he might 
express himself in any other language which he has mastered. 
And as with any other language, mastery requires much 
practice. 

Furthermore the student has two other important needs 
in the way of mathematical technique. These are 

1. How to integrate the differential equation. 

2. How to perform the calculations needed to obtain 
an answer correct within the practical requirements of the 
problem. 

4. The Fundamental Concept of Calculus. In chemical 
processes, and indeed in all processes involving changes, 
the first step is frequently to study the rate at which changes 
take place, and to analyse the causes of these changes. In the 
ordinary language of mathematics the expression " rate of 
change " is used for indicating the rapidity of change of 
one variable quantity as compared with another variable 
quantity. In discussing the rate of an engineering process 
we have to consider not only the quantity whose change we 
especially desire to analyse called the dependent variable, 
but also the independent variable, or quantity with which 
the dependent variable is compared. 

It is also important to distinguish clearly the two con- 
cepts " rate of change " and " average rate of change." For 
example, in chemical reactions the average velocity with 
which the reaction takes place is measured by dividing the 
quantity changed during a specified interval of time by that 
time. The rate of change at any instant is the derivative of 
the quantity with respect to the time, and is expressed as 
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-?, where dx represents the amount of the change of the 
dO 

quantity x during the time dO. 

While the expression " rate of change " is usually applied 
to changes involving time as the independent variable, there 
are a great many types of changes which do not explicitly do 
so, but which resemble time changes in character and be- 
haviour. It is convenient to speak of the rate of these 
changes also, the independent variable being necessarily 
specified in every case, and it will be found that the method 
of analysis used for studying time changes will apply equally 
well to the others.* 

5. Integrating the Differential Equation. The task of 
integrating the differential equation which occurs in a 
chemical problem may be simple or difficult, according to 
the nature of the problem. It is a happy circumstance that 
many of the most common and important problems are very 
easy in this respect, yielding differential equations of the 
" separable " type, or at worst are " linear "and therefore 
are already familiar to the student in proportion to his 
thoroughness in the calculus of the second year. A few 
harder cases will occur from time to time. 

A large proportion of these elementary types integrate 
in terms of natural logarithms. It would seem that Nature, 
on her physical and chemical side, has a peculiar fond- 
ness for the quantity e, the base of the natural system of 
logarithms, a fact which the mathematician may justify by 
remarking that e* is its own derivative. The student has 
usually had little or no practice in manipulating logarith- 
mic expressions and solving equations in which they occur. 
Such practice will inevitably be supplied in abundance if 
he works out the problems in the following chapters. The 
properties of logarithms as laid down in some good text- 

* For a fuller discussion of rates of change and of derivatives, consult 
Woods and Bailey's Elementary Calculus, Chapters I and II, or other 
standard textbook on calculus. 
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book on trigonometry or algebra may well be reviewed, but 
the explanations of the present text have been made very 
detailed and will probably be sufficient. 

6. Numerical Calculations. In regard to the im- 
portant matter of numerical calculation, there are three 
available methods in most problems, first, the use of tables 
including the slide rule, second, reduction of the problem to 
one in progression, and third, the use of graph paper. Usu- 
ally a happy combination of all these is best. As they are 
all taken up in detail in their proper place no further intro- 
ductory remarks are necessary, except to emphasize that, 
to the engineer as to the mathematician, the test of work- 
manship is in the soundness of his arithmetic. There is a 
real test of character in the answers a man will turn in. 
This test does not reside alone in their accuracy and freedom 
from blunders, essential qualities though these of course 
are. We may illustrate by citing two common cases of 
answers which we regard as vicious. One student has a 
question for which he is expected to return an answer cor- 
rect to three significant figures, the right answer being 75.4; 
he knows what is required of him and sets out with every 
fair and good intention to get it, but by, it may be, hasty 
work with a slide rule, he gets 73.2 and will then complain 
bitterly because he is given only a little credit for the prob- 
lem. Many a promising engineering job has failed to 
function for similar cause, but the excuse of this man usu- 
ally is that " chemical engineers don't need to know the 
answer better than five per cent." 

A second man, or more likely the same man next time, 
has a question where the experimental data by no stretch 
of interpretation can be correct to better than two signifi- 
cant figures. By zealous work with pencil or tables an 
answer such as 3.075398 is produced. Probably this result 
is more immoral than the former, since it has an appearance 
of accuracy which does not exist. 



8 DIFFERENTIAL EQUATIONS IN APPLIED CHEMISTRY 

The honest method of handling figures is simple and easily 
stated : 

1. Decide how many significant figures it is desirable to 
find. 

2. Decide whether the data warrant this or a smaller 
number. 

3. Make it a matter of religion to get exactly the number 
you have decided upon. 

7. Summary. The object of this introductory chapter 
lias been to point out the value of calculus as a tool for the 
chemist, and to emphasize that mastery of it requires exer- 
cise of the mental faculties in no immature way; first in 
analysing the problem and translating it into mathematical 
language, second in integrating or otherwise solving the 
resulting differential equation, and third in calculating the 
answer in a form suitable for practical use. The qualities 
of mind and the manipulative skill required for the three 
operations are quite distinct but all of value to the engineer. 
They are to be kept in view in all the work of this book. 



CHAPTER II 
PROCESSES OF THE FIRST ORDER 

8. Definition of a First Order Process. In many chem- 
ical reactions it is found that the rate of change of one 
substance into another is proportional to the first power of 
the concentration of the first substance. In these cases 
the velocity of the reaction at a given time depends only 
on the amount of the reacting substance at that time, pro- 
vided the other conditions of the reaction such as tem- 
perature, etc. remain constant. A reaction of this type 
is known to chemists as a reaction of the first order. And in 
other types of changes than chemical reactions it is frequently 
found that the rate of change of one variable as compared 
with another variable is proportional to the first power of 
the one variable alone. Such a change will therefore be 
spoken of as a process of the first order. 

To express a process of the first order in mathematical 
form, suppose that in a certain process involving the ma- 
terial x, the rate of change of x is proportional to the amount 
of x present. This would be a process of the first order. 

The rate of change of x (with respect to time) is -r-> as 

au 

pointed out in Chapter I. The process can therefore be 
expressed mathematically by the differential equation 

sr fa (1) 

where k is a " constant of proportionality " to be determined 
by some suitable means. 

The present chapter deals with processes of the first order. 
As a step toward a scientific understanding of natural proces- 

9 
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ses it is of the greatest importance to consider the question : 
why do so many of the commoner and simpler phenomena 
obey (at least approximately) an equation of the form (1)? 
The illustrative examples which follow below will suggest 
various answers to this question. Broadly speaking it may 
be said that equation (1) is a mathematical statement of 
the simplest case of that very general principle roughly 
known as the " Law of Mass Action." 

As an illustration, it appears to be true in the breaking 
down of radium atoms into smaller atoms, that for example 
out of one trillion radium atoms, the number which will 
explode per day is quite definite and fixed, independently of 
all external conditions whatever. It is evident that, if this 
be so, ten trillion atoms will yield ten times as many exploded 
particles per day as one trillion, or that, in general, the decay 
of radium follows the law (1). This also indicates that tho 
percentage decomposed per day is a constant, independent 
of the amount present, a fact that will be shown to apply 
mathematically, for all first order processes. 

Or again, take the case of a culture of bacteria with abun- 
dant food supply and no toxic conditions. If all other 
circumstances remain constant any group of a thousand 
bacteria will tend, by their natural law of growth, to in- 
crease in numbers just as fast as any other group of one 
thousand. The total rate of growth will, therefore, be pro- 
portional to the total number of bacteria present, in agree- 
ment with an equation of the form (1). 

Similarly, it is known in the theory of sound that the 
amount by which an organ pipe must be shortened to pro- 
duce a rise of a semi-tone in pitch is proportional to the 
length of the pipe, other things being equal. This is why, 
in a pipe organ, the tops of every set of pipes seem to lie 
on a curve of definite character. When equation (1) has 
been integrated and the character of the integral studied 
in detail, it will appear that the graphical representation of 
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a first order process of any kind is a curve which may be 
called " logarithmic " or " exponential " according to which 
of the variables we consider as independent. 

Example (A). It is known that the rate of inversion of 
cane sugar in dilute solution is proportional to the concen- 
tration of the unaltered sugar. (This was the result of 
experiment, of course.) If the concentration was 1/100 
when and was 1/250 when 6 = f> hours, express the 
concentration as a function of the time 6. 

Solution. Let the concentration of the unaltered sugar 
at time 6 be called c. Then 
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which is to be integrated under the assumption that the 
time is reckoned from the start of the experiment, that is, 
c = Co when = 0. By " separation of variables " (mul- 
tiplying by dd and dividing by c), the differential equation 
may be written: 

fir 

=kdO (3) 

c 

One way to solve this differential equation is to integrate 
both sides in the form 

In c = kO + const. (4) 

where " const." denotes the constant of integration, and 
where In c is the " natural logarithm " of c, that is, the 
logarithm to the base 2.71828 + denoted by e. (One reason 
why this base is properly called the natural base is that it 
enters inevitably and naturally in so many problems of 
physics and chemistry.) To determine the constant of 
integration we have the fact that c = Co when = 0, which 
gives by substituting in (4) 

In c = const. (5) 
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and by inserting this value of the constant in (4) we find 

In c = kO + In c (6) 

or, since the difference of the logarithms of two quantities 
is the same as the logarithm of their quotient, the equation 
may be written in the form 

ln--=/c0 (7) 

c 

(Equation (7) may also be obtained from (3) by taking the 
definite integral of both sides of (3) between proper limits, 
thus 

*-/'** (8) 

c Jo 

which gives (7) directly.) 

Now equation (7) holds true for all corresponding pairs of 
values of c and 0. Suppose that at a time 0i it was observed 
that c = d. Then 

In - = fc0! (9) 

Co 

Dividing (7) by (9) so as to eliminate ft, and partially clearing 
of fractions gives 



Co 01 Co \Co/ ^ 

where the last step follows because the logarithm, to any 
base, of a quantity raised to some power is the same as the 
index of the power times the logarithm of the quantity, 
an elementary and frequently used property of logarithms. 
Finally, since the logarithms of two expressions can be equal 
only when these expressions themselves are equal, 

c = co ( C Afl (11) 
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which is the desired relation between the concentration c 
and the time 6. Inserting the given data c = 0.01, ci = 
1/250, 0i = 5, 

c = 0.01 (0.4)1 (12) 

Remarks on the above solution. The equation (11) is a 
relation connecting the values of c at three different times. 
It may therefore be called the " three-point formula " for 
a first order process. 

The right member of (11) may be written in several other 
forms which are useful. For instance let 



M- 

= I - K = a 

W 



constant (13) 



equation (11) then becomes 

c = c r* (14) 

showing the analogy between a first order process and a 
geometric progression. 

It is to be noticed that, in working out the numerical 
solution of the above example, a literal solution was first 
obtained, the data afterwards substituted. This order of 
work is to be preferred when practicable for two reasons, 
first it yields a general result useful in all similar problems, 
and second it is usually quicker and more accurate. 

Example (B). Van t'Hoff found that 5.11 grams of di- 
bromsuccinic acid decomposed in hot water at the following 
rate: 

When = 10 20 30 40 50 60 minutes 
Acid = 5.11 3.77 2.74 2.02 1.48 1.08 0.80 grams 

Set up an equation of the form (2) and calculate k for each 
of the times given. Is k sufficiently constant to indicate 
a reaction of the first order? 
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Solution. Letting c stand for the number of grams acid 
remaining at time 6, integrating as in Ex. (A), and solving 
equation (7) for fc, 

k = (Inc - In Co) -*- (15) 

from which fc may be calculated by a table of natural loga- 
rithms. To use common logarithms advantage is taken of 
the relation 

In c = 2.303 log c (16) 

where log denotes the logarithm to base 10. Thus 

A: = (log c - log co) (2.303) -s- 9 (17) 

Using four-place logarithms and arranging the calculation 
in tabular form we find: 

log co = log 5. 11 = 0.7084 



10 
20 
30 
40 
50 
60 

The calculations for the last three observations arc left to 
the student. 

Example (C)r A tank contains 100 gallons of brine, the 
amount of dissolved salt being 50 pounds. Fresh water is 
allowed to run into the tank at the rate of three gallons per 
minute, the brine overflowing at the same rate. The solution 
is kept practically uniform by agitation. How much salt 
will there be in the tank at the end of one hour? 

Solution. The simplest way to set up a differential equa- 
tion for problems on rate of flow is to state the obvious 
relation 

increase = inflow minus output (18) 



logc 


- log (c/c ) 


- In(c/c ) 


- k 


0.5763 


0.1321 


0.3041 


0.03041 


0.4378 


0.2707 


0.6232 


0.03116 


0.3054 


0.4031 


0.9281 


0.03094 
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Considering an interval of time dd the relation (18) may be 
stated more at length thus: the increase dc which takes 
place in the amount of salt in the tank during time dd is the 
difference between the amount of salt flowing in during 
this time and the amount flowing out. This fact is axiom- 
atic provided no chemical reaction or other cause creates 
or destroys salt inside the tank. It is equally evident that 
inflow = (rate of flow of solution) (concentration) (time) (19) 
with a similar relation for output, and also 

concentration = total salt -=- volume (20) 

Calling V the volume of brine, c the total salt at time 0, 
it is clear that in this problem the inflow of salt is zero since 
the incoming liquid is fresh. The concentration of the out- 
going solution is c/ V arid the rate of outflow of solution is 
3 gal./min. Hence 

output of salt = 3y dO (21) 

meaning that in an interval of time dO short enough so that 
during this interval the concentration does not sensibly 
change the number of pounds of salt which leaves the tank 
is 3cdO/V or 0.03 cdO since the volume of brine is con- 
stant at 100 gals. Substituting in (18) 

dc = -0.03cd0 (22) 

Integrating in the form (7) gives 

ln^= -0.030 (23) 

ou 

To solve this equation for c, use the important property 
of logarithms 

e ln = q (24) 
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where q is any quantity whatever. (This is merely one 
form of the definition of a logarithm.) Therefore from (23) 

c = 50 e - 0.030 (25) 

Putting 6 = GO min. and using a table of exponentials gives 
c = 8.265 Ibs. (Another good way is to write (23) in the 
form 

2.3 (log c - log 50) = - 0.03 (60) (26) 

then solve for log c and use a table of common logarithms.) 
Example (D). The air in a recently used class-room 
30 X 30 X 12' tested 0.12% by volume of carbon dioxide. 
How many cubic feet of frosh air must be admitted per 
minute in order that ten minutes later it shall contain not 
over 0.06%? (Fresh air contains 0.04% CO L >.) Assume 
immediate mixing of frosh with stale air. 

Solution. Set up equation (18). Adopt the notation 

y = concentration of CO 2 by vol. at time 

a = cu. ft. /min. of entering air. 

V = vol. of room. 

?/o = concentration CO 2 at start. 

g = concentration CO 2 in fresh air. 

It is then evident by (19) that, in time d6, 

inflow = ag dd (27) 

and, since stale air leaves by leakage through doors and 
windows at the same rate at which fresh air enters, 

output = ay dd 28) 

^ 
The total CO 2 at time is Vy, by a relation like (20). Hence 

the increase in time dO is d (Vy) or Vdy, since V is constant. 
Filling in equation (18) 

Vdy = agdS - ay d9 (29) 
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The easiest way to solve this differential equation, and at 
the same time to see that the process is of the first order, 
is to introduce a new variable; x defined by 

x = y- (30) 

whence dx = dy because g is a constant and drops out on 
differentiation and (29) may be written 

dx a 



which is of the form (1) and shows the value of k to be 



> that is 



a = - kV (32) 



making it clear that the unknown a is found when we find 
k. Proceeding as in Kxample (A), equation (8), 



dx/x = k C 9 dd (33) 

/() 



and solving for k as in (15) 

k = (In X/XQ) + (34) 

or by (32) 

a = - V (In ar/xo) -5- (35) 

which by (30) is the same as 

a=-^ln^^- (36) 

6 2/o - g 

We now insert the data of this example, V = 10,800, 0=10, 
y = 0.0006, 7/0 = 0.0012, g = 0.0004 and have 

10,800 . 0.0006 - 0.0004 inen . t ,. 

a - ~ -10- ln 0.0012 - 0.0004 = ~ 108 ln 1/4 
= + 1080 In 4 = 1080 (1.386) = 1500 cu. ft. apprax. 
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Remarks on the above solution. The process here em- 
ployed of introducing an auxiliary variable to reduce the 
differential equation to the form (1) is very important. Its 
significance will appear more fully on considering the two 
important properties of first order processes next to be given : 

1. The geometrical progression property. 
2. The semi-logarithmic property. 

9. The Geometrical Progression Property. It has already 
appeared in equation (14) that a formal resemblance* exists 
between a first order process and a geometric progression. 
This becomes a real resemblance when the exponent 6 is an 
integer. For any term of a geometric, progression is equal 
to the first term multiplied by an integral power of the 
ratio r. The analogy may be more* fully stated thus: 

If a quantity varies ace;ording to a differential equation 
of the form (1), and if in some definite interval of time the 
quantity is known to e^hange by a e*ertain fraertion e>f its 
value, it will change by the same fraction during any other 
equal interval of time. 

This fact is most freejuently used by stating the time in 
which a variable, obevying equation (1), falls to half its 
value. For instance the rate of decay of the radioactive 
substance thorium X is such that the activity falls to half 
its value in four days. Or the length of an organ pipe is 
halved for every octave of rise in pitch, other things being 
equal. But there is no special virtue in the fraction one- 
half except, simplicity, and the principle is of great use in 
calculation. 

Example (E). Using the result of Example (C), find how 
much salt remains in the tank after two hours. 

Solution. Instead of putting 6 = 120 and solving (25) or 

83 

{26) we note that the salt fell to ^^r of its value in one hour. 

oUU 
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Therefore in two hours its value will be 



which is readily computed by a slide rule. 

Example (F). With everything as in Example (D), 
find what would be the concentration of C() 2 if the same 
current of fresh air could be continued for a second period of 
ton minutes. 

Solution. It is important to notice first that it is the 
variable x and not ?/ which satisfies the equation of form 
(1). Now x at the start was ?/ - g that is 0.0012 - 0.0004 
or 0.0008. After ten minutes x has fallen to ().()()()() - 0.0004 
or 0.0002. Thus x lias fallen to \ its value in 10 mill. By 
the G. P. property x will in another ten minutes fall to }- 
of its value 0.0002 or 0.00005, and ij will become 0.00045 
or 0.045 %. 

Remark on the last two solutions. In practice it often 
happens that problems on first order processes can be solved 
mentally or with slight labor using the G. P. property, which 
would be otherwise lengthy . This will usually be true when 
the time intervals which occur are multiples of some common 
divisor. 

10. The Semi-logarithmic Property. One of the most 
convenient ways of studying a process is to devise some 
method of plotting on graph paper so that the graph shall 
be (approximately) a straight line. If the process is of the 
first order there are two excellent methods of doing this. 

First Method. By inspection of equation (17) it appears 
that if we plot log c as ordinate against 6 as abscissa on 
ordinary graph paper, the result will be a straight line, 
providing always that c and 6 are connected by the differ- 
ential equation (1). The data of Example (B) may therefore 
be easily tested to see if the process is of the first order by mak- 
ing such a plot, without actually calculating the constant k. 
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The points obtained will lie on a straight line, due allowance 
made for experimental error and for possible actual depar- 
ture of the phenomena from a true first order process. The 
verification is left to the student. When making the plot, 
the following suggestions should be borne in mind: 

1. Choose a piece of graph paper 6" X 6" or larger, 
ruled decimally whether in centimeters or inches is a 
question of taste. 

2. Make a table of the values to be plotted. In the data 
of Example (B), the first column 6 is abscissa against the 
second column log c as ordinate. 

3. Make a horizontal and a vertical scale 1 such that the 
largest and the smallest values to be plotted fall near the 
respective ends of the scale. The full extent of the scale 
will thus be utilized. 

4. In marking points, estimate tenths of a scale division 
by the eye, and neglect values smaller than a tenth. Points 
are best indicated by small finely drawn + signs. 

Second Method. It is often convenient to plot the same 
graph by using semi-log paper. We plot c as ordinate 
against as abscissa, using for c the scale already printed on 
the graph, with suitable introduction of the decimal point. 
Paper may be obtained with either one, two, three or four 
log scales to the page. The graph is continued across the 
page as many times as necessary, each journey across the 
page being numbered, and a corresponding scale indicated 
on the margin. It is evident that these successive passages 
over the page will indicate a set of parallel lines. It is thus 
possible to graph on a single page a wide range of the vari- 
ables more conveniently than with ordinary paper, though 
sometimes with less accuracy. 

The method of plotting on semi-log paper is illustrated 
in Fig. 1, which is a graph of Example (A), carried from 
6 = to B = 18 hours. During this time the concentration 
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of the unaltered sugar falls from 0.01 to 0.00037. The hori- 
zontal scale is divided into six equal intervals, each corre- 
sponding to one hour, so that three journeys across the page in 
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Semi-log Plot of Ex. (A) 
FlO. 1. 

the horizontal direction are required to cover eighteen hours. 
The vertical scale is logarithmic. On the left of the diagram 
the decimal point has been so placed that we may begin at 
the top with the concentration 0.01, corresponding to 
6 = 0, and end at the bottom with the concentration 0.001, 
corresponding to approximately 6 = 12.6. On the right 
of the diagram the decimal point has been moved one place 
to the left, giving 0.001 at the top and 0.0001 at the bottom. 



22 DIFFERENTIAL EQUATIONS IN APPLIED CHEMISTRY 

The graph consists of four straight lines numbered in order. 
The first covers the time from to 6 hours, the second from 
6 to 12 hours, while the third and fourth, which together 
make one horizontal journey across the figure, cover the 
time from 12 to 18 hours. On the other hand, the first, 
second, and third pieces of the graph together make one 
journey from top to bottom, and cover the change in con- 
centration from 0.01 to 0.001. The fourth piece is part of 
a second journey in the vertical direction, beginning at 0.001 
and ending at 0.00037. The plot was made by first noting 
the two observations given in the text, namely c = 0.01 
when = 0, and c = 0.004 when 9 = 5. These are in- 
dicated by two small circles about the points in question. 
A straight line is drawn through those two points and con- 
tinued with constant direction as described above. 

It should be carefully noticed that in extrapolating, that 
is predicting the course of a process over an interval of 
time extending far beyond the range of the observed values, 
three sources of error exist : 

1. Errors in drawing (or calculating). These arc likely 
to be magnified by extrapolation. 

2. A slight error in the slope of the straight line, due to 
unavoidable experimental errors in the data themselves, 
may produce a relatively large error in the predicted values. 

3. The actual process almost invariably departs more 
or less from the assumed law. Such departure may be 
gradual and progressive, or may be very sudden, according 
to the nature of the phenomenon. 

PROBLEMS 

1. With the data of Ex. (A), find the value of c after 13 hours. 
(Slide rule accuracy is sufficient.) 

2. Show that the answer to Ex. (A) may be written 

c = 0.01 e -- 1830 

When would this form of answer be convenient? Verify the answer 
to Prob. 1 by use of this formula. 
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3. Make the table for = 40, 50, 60, in Ex. (H), using four-place 
logs as in the text. When four-place logs are used, to how many sig- 
nificant figures will the calculations be accurate? 

4. Find the average of the six values of k in Ex. (B). What is the 
greatest per cent departure from this mean value? What is your con- 
clusion in regard to the phenomena? Can any rule be given to 
govern your decision as to whether the process is of the first order? 
Is this average value of k the "most probable value?" 

5. With the data of Ex. (A), let z denote the amount of invert sugar 
at time 0. Express z as a function of 0. 

6. In a first order reaction docs the value of k depend upon the 
amount of substance present at the start? Does it depend on the 
unit of mass employed? Does it depend on the unit of time em- 
ployed? 

7. With the notation of Prob. 5, suppose we were able to observe 
Co, zi, 2, 0i and 2 . Express k in terms of these quantities. 

8. With everything as in Ex. (C), after how many hours will the 
solution contain less than half of one per cent salt by weight? (Slide 
rule accuracy.) 

9. A student working Ex. (C) made the mistake of using equation 
(11), taking 9i = 1, and ci = 48.5. If his calculation were otherwise 
correct, what was the per cent error in his answer? (Use four place 
logs.) Can you state why his logic was bad? 

10. Blariksma, when transforming 49.30 gms. acetochloranilid 
into para-chloracctanilid, observed that 13.70 gms. of the para com- 
pound were formed during the first hour, and 9.85 gms. during the 
second hour. What evidence do these data give that the process is 
of the first order? (Use G. P. property.) 

11. With the data of Ex. 10, how much acetochloranilid would we 
expect to find unchanged after 8 hours? (Use the C. P. property.) 
lie found actually 4.8 gms. How do you explain the discrepancy? 

12. Make a graph of Prob. 10 by plotting the log of the acetochlor- 
anilid against the time on ordinary cross-section paper. Read from 
your graph the number of hours needed for the acetochloranilid to fall 
to 0.1 gm. How accurately can the time be read? How well does the 
value of the acetochloranilid read from the graph when = 8 agree 
with the calculated value from Prob. 11? Compare your reading for 
0.1 gm. with the rest of the class. 

13. Make a graph of Prob. 10 by plotting acetochloranilid against 
time on semi-log paper. How well can you read the answer to Prob. 
11? Read also the time needed for the acetochloranilid to fall to 
0.1 gm. 
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14. Using the data of Ex. (B), plot log c against as described in 
the text. Draw a straight line which shall pass approximately through 
the seven points, and continue it far enough to read the value of log c 
when 6 = 100. Repeat using semi-log paper. 

The following problems arc to be solved by first setting 
up a differential equation which, in each case, may be thrown 
into the form (1) or (2). The solution may then be com- 
pleted by using one of the integrated forms, like (7), (11), 
(17), or (25), or by using the G. P. principle, or by making 
a graph. The student should not limit himself to a single 
method, but should select in each case what appears the most 
convenient. It is well to check the answer by using two 
different methods. Slide rule accuracy will be sufficient. 
Ordinary data such as specific gravities may be found from 
a handbook when needed. 

16. If the activity A of a radioactive deposit is proportional to its 
rate of diminution, and is found to decrease to -J its value in 3 hours, 
what is its value after 9 hours? After 13 hours? 

16. According to Newton's law, the rate at which a substance cools 
in air is proportional to the difference between the temperature of the 
substance and that of the air. If the temperature of the air is kept 
at 20 C., and the substance cools from 100 to 60 in 20 minutes, when 
will its temperature become 30? When 21? 

17. A porous material dries in the open air at a rate approximately 
proportional to its moisture content. (Can you see a reason why this 
should be so?) If a sheet hung in the wind loses half its free moist- 
ure the first hour, when will it have lost 99 per cent, weather conditions 
remaining the same? 

18. If in the culture of yeast the amount of active ferment doubles 
in 3 hours, how much may be expected at the end of 15 hours at the 
same rate of growth? At the end of 20 hours? 

19. Hydrogen flows through one tube into a litre flask filled with 
air and escapes through another. If the current is so slow that the 
mixture may be considered uniform, what percentage of hydrogen will 
the flask contain after 3 litres of gas have passed through? How much 
after 10 litres? 

20. Repeat Prob. 19 using oxygen in place of hydrogen. Assume 
that air contains 21 per cent of oxygen. 
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21. If the average person takes 18 breaths per min., exhaling each 
time 100 cu. in. containing 4% carbon dioxide, find the per cent CO* 
in the air of a class room \ hour after a class of 50 enters. Suppose the 
air fresh at the start, and the ventilators admit 100 cu. ft. fresh air 
per min. (Fresh air contains 0.04% CO2 by volume.) Take volume 
of room = 10,000 cu. ft. 

22. A factory room 200 X 45 X 12' receives through the ventilators 
10,000 cu. ft. fresh air per min., which contains 0.04% of C02 by vol. 
The help enters at 7:00 A.M. A half hour later the CO2 content has 
risen to 0.12%. What value is to be anticipated at noon? How much 
fresh air must be admitted per min. to insure that the COa shall never 
rise above 0.10%? 

In such a problem will the ventilation be better or worse if the fresh 
air and the foul air do not mix rapidly? If we assume no mixing at 
all, what will be the differential equation? 

23. A lot of leather is washed free from salt which it contains, by 
suspending it in water in a drum, into which fresh water is flowing at 
the rate of two gallons per minute. The leather is very thoroughly 
agitated by rotating the drum and the solution of salt thus formed 
runs to waste through an overflow pipe as fast as fresh water is admitted. 
The leather is so thin and porous that the diffusion of the salt into the 
surrounding solution may be considered as instantaneous in comparison 
with the rate at which the water is being changed. 

The drum contains 500 gals, of water. Each charge contains 400 
Ibs. of leather on a dry basis; this when taken from the drum and sent 
to the dry loft contains 600 Ibs. of solution in its pores. If at the start 
of the operation each charge contains 100 Ibs. of salt, while the dry 
product must not exceed 0.4% how long must the washing be con- 
tinued? 

24. The population of New England increased from 26,000 in 1640 
to 1,000,000 in 1790, during which interval there was very little immi- 
gration. Allowing 30 years to a generation, estimate the average 
number of children per family, no allowance being made for infant 
mortality nor unmarried persons. 

25. The population of the country is now 100 millions. The nat- 
ural growth of population would cause the number of inhabitants to 
double in fifty years if there were no immigration, which, however, 
brings a constant influx of one million per year. Assuming the 
natural increase among the immigrant population to be no higher 
than among the native element, in how many years will the popula- 
tion double? 

26. If a quantity increases according to a first order process, what 
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must be the value of k in order that the rate of growth may be equiva- 
lent to 6% interest compounded at the end of every time-unit? 

27. Prove that continuous compound interest at 5% is very nearly 
equivalent to yearly compound interest at 5.1%. 

28. If a man can earn 5 dollars a day over expenses, and keep his 
savings continuously invested at 6% compound interest, how long will 
it take him to save $25,000? 

29. Consider a vertical column of air, and assume that the pres- 
sure at any level is due to the weight of air above. Show that p = p Q e-i* 
gives the pressure at any height h, if Boyle's law be used. Determine 
k from the known properties of air, assuming the temperature of the 
whole atmosphere to be C. 

30. A first order process gave a yield of one gm. during the first 
hour and 0.7 gm. during the second hour. What was the total yield? 

31. In every gram of radium, 13 X 10 10 atoms change per sec. 
into the emanation. Find how long for J the radium to be changed. 
Given that 1 c.c. of hydrogen contains 2.7 X 10 19 molecules, and the 
atomic weight of radium is 226. 

Ans. about 1900 years. 

32. The element thorium, which is not itself radioactive, changes 
into active thorium X at a substantially constant rate. The active 
thorium X may be separated from its parent thorium by precipitation 
with ammonium hydroxide, the active part remaining in the filtrate. 
After the filtration it is found on the one hand that the activity of the 
filtrate falls to J its value in 4 days, and, on the other hand the precipi- 
tate recovers activity on standing. Find how long before the precipitate 
will recover \ the former activity. 

33. With the data of Prob. 32, what fraction of its original activity 
will the inactive thorium recover in 1 hour? In 5 days? 

34. Some radium chloride in solution, of normal activity 1000, 
was deprived of its emanation by a rapid current of air. It was then 
allowed to stand air-tight for 105 mins., when its activity was found to 
be 13.1. Find the radioactive constant k for the decay of the emanation. 



CHAPTER III 
PROCESSES OF THE SECOND ORDER 

11. Definition of Second Order Processes. In the usual 
language of chemists, a reaction is said to be of the second 
order when the rate of the reaction is proportional to the 
product of the concentrations of two reacting substances, or, 
as a special case, to the square of the concentration of one 
reacting substance. To say the same thing in mathematical 
symbols, let c and c f be the concentrations of the two react- 
ing substances, and x the concentration of ono of the prod- 
ucts. If the reaction proceeds in agreement with the 
differential equation 

fe - kcc ' < 37 > 

it is said to be of the second order. 

If a and b denote the amounts of the two substances 
present at the start, and if our unit of measure is the gram- 
equivalent of each substance, it follows that 

c = a x y and c' = 6 x (38) 

by which (37) takes the form 

f?T 

jg = fc (a - x) (b - x) (39) 

This differential equation may be described in mathe- 
matical language by saying that the derivative of the de- 
pendent variable with respect to the independent variable is 
equal to a polynomial of the second degree in the dependent 
variable. 

27 
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Many processes other than chemical reactions proceed in 
agreement with a differential equation like (39). By anal- 
ogy with the language of the chemist, these will be called 
processes of the second order. 

Remark. It must not be said, however, that the differ- 
ential equation itself is of the second order, for to the math- 
ematician this would mean the presence of the second 

d-y 
derivative ~ That a word shall be used in different 

senses by different people and at different times appears a 
law of nature difficult to alter, a fact, at least, which is true 
of the word order even within the domain of mathematics. 
No confusion is caused, because a mathematician is ex- 
ceedingly careful about his language. In this book, in- 
tended not so much for mathematicians as for chemists 
who wish to absorb mathematical ideas, it appears better 
to employ a term in a way familiar already. The two 
expressions " process of the second order " and " differential 
equation of the second order " will cause no confusion, par- 
ticularly since there will be little occasion for use of the 
latter, more mathematical, term. 

Differential equation's like (39) occur very frequently in 
that type of chemical reaction known as metathesis. 

Example ((7). Reicher found that 0.5038 gram mols of 
sodium hydroxide acted on 0.3114 gram mols of ethyl acetate 
at the following rate: 

When 6 = 393 G09 1010 1205 

alkali = 0.5638 0.4806 0.4467 0.4113 0.3879 

ester = 0.3114 0.2342 0.1943 0.1589 0.1354 

Show that the reaction is of the second order. 

Solution. The chemical equation is as follows: 

NaOH + C 2 H 5 Ac = C 2 H B OH + NaAc 

On the assumption that the reaction is substantially non- 
reversible, the rate of formation of the alcohol would be 
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proportional to both the concentration of the hydroxide 
and that of the ester. If a and 6 denote respectively gin. 
mols of alkali and of ester at the start, and x the mols of 
alcohol at time 6, then a x and b x will be the mols 
alkali and ester left to react at time 0. Hence x will satisfy 
(39) provided this conception of the nature of the reaction 
is true. 

To try it out, proceed as in Example (B). Integrate 
(30) between appropriate limits, 

r* * = r** kde ( 4o) 

Jo (a - x) (b - x) Jo v ' 

yielding the result 
2.302(>riog (a - x) - log (b - x) - log-^1 = (a - 6) kO (41) 

By making a table after the style of Example (B) it will be 
seen if k is constant within the limit of experimental error. 
The trial is left to the student. Four-place tables should be 
used. 

Example (//). In a laboratory experiment, spent iron 
oxide from a gas works, containing 52% of sulphur, was 
extracted with boiling benzol. The following data were 
obtained: 

Weight of spent oxide = 25 gms. 
Weight of benzol =100 gms. 

Time in 

minutes 10 20 30 40 50 60 70 80 90 

Concen- 
tration 2.15 3.4 4.32 5.10 5.73 6.32 G.7S 7.04 7.42 

where the concentration is expressed as grams of sulphur 
dissolved per 100 grams of benzol. 

The solubility of sulphur in benzol at the boiling point is 
11.7 grams per 100 grams of benzol. 
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It was desired to fit an equation to the data which should 
represent the facts within the precision of the experiment, 
about 5 per cent. 

Solution. It was first assumed that the rate of solution 
of the sulphur was proportional to the amount of sulphur 
left undissolvcd, which would imply a process of the first 
order. To test this, the logarithm of the number of gins, 
of undissolved sulphur was plotted against time. It was 
found that the best graph through the points was a curve, 
and not a straight line* as it would have been had the process 
been of the first order. 

It was then assumed that the rate of solution of the 
sulphur was proportional both to the amount of undissolved 
sulphur and also to the difference between the concentration 
of the solution and that of a saturated solution at the same 
temperature. This is actually a more reasonable suppo- 
sition than the former one, because no sulphur would dis- 
solve when the solution became saturated, even though 
there remained some undissolved. If the improved con- 
ception of the process is true it will be of the second order 
and the differential eq nation will be 

g = - kw (S - c) (42) 

where w is the number of grams of undissolved sulphur at 
time 9, while 8 and c are the grams of sulphur per 100 grams 
boiling benzol in a saturated solution and in the actual 
solution at time 0, respectively. 

This differential equation contains three variables, w, c, 6. 
But w arid c are connected in virtue of the fact that 

Total sulphur present = c + w = 25 X 0.52 = 13 grams (43) 
Hence w = 13 c and dw = dc giving 

= + * (13 - c) (11.7 - c) (44) 
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an equation of the form (39). By integrating as in (41) 
2.3026^ (13 - c) - log (11.7 - c) - log ^~J 



= (13 - 11.7) k6 



(45) 



If this equation is a good representation of the process, a 
plot of the quantity log [(13 c) /(1 1.7 c)] against 8, 
using the observed data, should yield a series of points nearly 
on a straight line. 
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Remarks on the above solution. It is interesting to note 
(Fig. 2), that what appears to be the best straight line 

through the points does not pass through the pointfo, logd> 
a fact probably due to a considerable error in observing 
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the start of the experiment. The actual start seems to 
differ by about three minutes from that indicated by the 
first observation. 

The two assumptions made as to the rate of solution of 
the sulphur appear very reasonable when we try to under- 
stand the mechanism of diffusion. On the surface of the 
undissolved sulphur is a thin film of saturated solution. 
The forces which carry this film off into the neighboring 
dilute solution are analogous to a difference of potential. 
That the rate of solution should also be proportional to the 
amount undissolved is equally natural, for when the soluble 
material is contained in a mass of inert substance the sur- 
face exposed is substantially proportional to the amount 
present. 

12. Finding the Best Line through a Set of Points. 
In Example (H) it was seen that the plot of log w against 6 
gave on the whole a curve and not a straight line. The plot is 
shown on Fig. 2, where the scale on the right of the page 
corresponds to log w and the observed points are indicated 
by the mark X . 

13 c 

The scale on the left of Fig. 2 corresponds to log -r-^ 

1 1. / c 

and the observations are indicated by the mark + These 
points lie on a straight line, within the limits of experimental 
error. 

This procedure of devising some method of plotting by 
which a series of points are made to lie nearly on a straight 
line is very useful, and is the method usually employed in 
fitting an equation to a set of data. It now becomes im- 
portant to draw the best straight line through these points. 
A rough method is merely to draw the line with a ruler, 
estimating the position by the eye so as to leave about as 
many of the observed points on one side of the line as on 
the other. 

A slightly more accurate method, and the one used in 
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plotting Fig. 2, is called the Method of Averages. It is shown 
in works on the Theory of Probability and Method of Least 
Squares that the most probable straight line through a set 
of points (under certain assumptions as to the nature of 
the experimental errors) passes through the point whose 
abscissa is the average of the abscissas of all the points of 
the set, and whose ordinate is the average of all the ordi- 
nates. That is, to find this Average Point, we add together 
the abscissas (in Example (H) these will be the values of 8) 
and divide their sum by the total number of points in the 

group. Similarly we add the ordinates, (in Example (H) these 

13 c 

will be the values of log 7 ), and divide the sum by the 

JL JL / C 

total number of points in the group. The point whose 
abscissa and ordinate are the numbers thus found will be 
called the Average Point of the Group. 

Since two points are needed to determine a straight line, 
the whole observed set of points is divided into two groups, 
and the average point of each group is found. In Example (H) 
are ten observations. The first five and the last five, respec- 
tively, were taken, and the average points found. These 
two points are indicated at AI and A 2 , Fig. 2, by small 
circles about the points. Finally the straight line was 
drawn through these two points by a ruler. 

The method of averages is usually quite as rigorous as 
the data warrant. The equation of the straight line through 
the two average points may be calculated by analytical 
geometry, (Woods and Bailey, Elementary Calculus, Art. 14). 
Hence the value of the constant k and the corrected reading 
for the time at the start of the experiment may be found 
with greater accuracy than would be possible from only 
two observations. For fuller discussion, works on the 
method of least squares may be consulted. For other 
graphic methods, and many examples, see Lipka's Graphical 
and Mechanical Computation, Chap. VI. 
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PROBLEMS 

1. Make a table for Ex. (G) after the manner of Ex. (B). Hence 
find the average value of k for Reicher's experiment. 

2. In a second experiment Reicher began with 0.3910 gm. mols 
alkali and 0.6593 gin. mols ester, and found after 342 minutes 0.2885 
of alkali and 0.5568 of ester. Calculate the value of k. 

3. In Prob. 2 name two conditions essential to the experiment if k 
is to be the same in both runs. 

4. Derive equation (41) from (40). (Lipka's table of integrals may 
be used.) 

6. Solve equation (41) for x. Hence find the amounts of sodium 
hydroxide and of ethyl acetate to be expected when 6 = 684 in the 
second experiment of Reicher (Prob. 2). 

6. Make a graph for the data of Prob. (2) by plotting log : 

U X 

against on ordinary cross section paper. Hence verify the result 
of Prob. 5. 

7. Make a graph for the data of Prob. (2) by plotting = against 

o x 

8 on semi-log paper. Hence verify the result of Prob. (5). 

8. Arrange equation (41) so as to show an analogy with geometric 
progression. State in words the G. P. property of a process of the 
second order. Use this principle to verify the result of Prob. 5. 

9. Use your graph of Prob. 6 to determine the value of k, and com- 
pare the result with Prob. 2. 

10. Suppose Reicher had used equivalent mols of sodium hydroxide 
and ethyl acetate, for instance 0.500 mols of each to start, other con- 
ditions remaining as before, how long would have been required for 
half the ester to be decomposed? 

11. In the special case of a second order process where a = 6, how 
may the data be plotted so as to give a straight line? Using both the 
data and the result of Prob. 10, make such a plot, and read off the 
value of the time when the ester shall have fallen to i its original 
value. 

12. State in words the property possessed by the special case treated 
in Probs. 10 and 11 which enabled you to plot a straight line. Hence 
verify the result of Prob. 11 by a simple calculation. 

13. Compare the course of a chemical reaction of the second order 
with one of the first by plotting the amount of product against time. 
Use the same axes for both reactions, and assume initial amount of 
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reacting substances in each case to be one gram mol. and each re- 
action to be half completed in one hour. 

14. If data are available showing the yield of a process for a number 
of values of the time, and if the initial amount of all the reacting sub- 
stances is known to be one gram mol., what do you regard as the 
easiest way to determine whether the reaction is of the first order, or 
of the second order, or neither? Test the following: 



(a) 


When 6 


= 


1.0 


1 


.7 


2.5 


5.0 


10. 





<x 


hours 




X 


= 


0.20 


0. 


,29 


0.40 


0.50 


0. 


70 


1.00 


mols 


(b) 


Whenfl 


= 


1.0 


1. 


4 


3.3 


5.0 


10. 





oc 


hours 




X 


= 


0.20 


0. 


22 


0.25 


0.33 


0. 


68 


1.00 


mols 



(c) When = 05 8 12 20 oc hours 

x = 0.50 0.68 0.82 0.95 1.00 mols 

15. Make a graph of Ex. (G) in the style of Fig. 2. Calculate the 
value of k by the method of averages. To how many decimal places 
is it well to compute fc? Compare your answer with the result of 
Prob. 1. 

16. A process is known to satisfy a differential equation of the form 
(39). The value of the constant b is known to be 9.82 mols. The 
constant a cannot be directly measured. The following observations 
of x were made: 

When = 20 40 * hours 

x = 5 7.76 9.82 rnols. 

Find the constants a and A'. 

17. A mass of inert material containing 5 Ibs. of salt in its pores 
is agitated with 10 gals, of water. In 5 minutes 2 Ibs. of salt have 
dissolved. When will the salt be 99% dissolved? 

18. A mass of insoluble material contains 30 Ibs. of salt in its pores. 
The mass is agitated with 20 gals, of water for one hour, when \ the 
salt is found to be dissolved. How much would have dissolved in the 
same time if we had used double the amount of water? (Assume the 
mechanism of solution to be like the extraction of sulphur in Ex. (H). 
A saturated solution of salt holds 3 Ibs. per gal.) 

19. A porous material dries in a confined space as if the moisture 
it contains were dissolved out by the air. A quantity of material con- 
taining 10 pounds of moisture was placed in a closed storeroom of vol- 
ume 2000 cubic feet. The air at the beginning had a humidity of 25%, 
(that is, } saturation). At the temperature of the storeroom saturated 
air holds 0.2 Ib. moisture per Ib. of dry air. If the material lost half 
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its moisture the first day, estimate its condition at the end of the 

second day. 

(A cubic yard of air weighs approximately two pounds.) 

20. How long would be needed for the substance of the preceding. 

problem to lose 90% of its moisture if the humidity of the air is kept 

at 25% by ventilation? 



CHAPTER IV 
SIMULTANEOUS PROCESSES 

13. Side by Side Reactions. There are many proces- 
ses both in chemistry and in engineering which may be most 
simply investigated by regarding them as the result of two 
or more simpler processes going on simultaneously. The 
point of view will be made clear by taking up a number of 
special cases. 

Side Reactions 

The expression " side reactions v is an abbreviation for 
the longer expression " side by side reactions." 

Example (I). A certain radioactive substance yields 
two quite distinct emanations, each at a different rate. 
It is known that the rate of formation of each emanation is 
proportional to the amount of the parent substance present. 
Each emanation is formed directly from the parent. Thus 
the process consists of two simple first order reactions going 
on side by side. 

Call c = amount of substance present at time 6. 

x = amount of first emanation at time 0. 

y = amount of second emanation at time 6. 
fci = velocity constant of first reaction. 
fc 2 = velocity constant of second reaction. 

The differential equations for the separate reactions may 
then be written 

|. fee and g-J* (46) 

37 
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These equations contain the four variables x, y t c t and 8. 
But x, y, and c are related in virtue of the fact that 

x + y + c = Co, the original amount of the parent. (47) 
The following data were obtained: 

When 6=0 3 days 

c = 10 5 units 

x = 4 " 

y= 1 " 

It is required to find the values of k and k z and to express 
x, y, and c as functions of the time 0. 

Solution. An easy way to begin is to divide the two 
differential equations one by the other, so as to eliminate 
dO. This gives 

? = T (48) 

dx ki ^ ' 

fo> 

By separating variables this is the same as dy = j^dx. 

K\ 

Integrating, and noticing that y is zero when x is zero, it 
follows that 

-K- 1-1 < 49) 

that is, the ratio of the two products of the process is constant 
and equal to the ratio of the velocity constants, a characteristic 
property of " side reactions " when the two or more simple 
reactions are of the same type, known as " Wegscheider's 
principle." This principle is very important, since it helps 
in practice to distinguish side reactions from reversible or 
consecutive reactions, to be considered later. Now since 
x = 4 when y = 1 we must have 

fci = 4 fe (50) 
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An easy next step is to add the two differential equations 
(46), giving 

^ - ft + to c (51) 



but by differentiating (47) 

d (x + y) = - dc (52) 

whence (51) becomes 

^ = - (fa + fe) c (53) 

which is of the form (1). In other words the parent sub- 
stance is used up as if the process were a simple first order 
reaction with velocity constant equal to the sum of the 
actual velocity constants. Therefore an equation of the 
form (15) may be at once written down, namely 

0=.-ln^-8-3 = |ln2 (54) 

4 1 

By (50) it follows that ki and fe are respectively = and = 

o o 

their sum. Hence fa = I In 2 = = (0.693) = 0.185 and 

O o lO 

fc 2 = 0.046. 

The easiest way to express c as a function of is to use the 
" three-point " formula (11), 

/\\e 

c = 10(5)8 (55) 

W 

Finally we see from (47) that x + y = 10 c and from (49) 
that x and y are respectively 4/5 and 1/5 of their sum. 
Hence 



and 



j/ = i[lO - 100)1] (56) 
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Example (/). In preparing dinitrobenzene from mono- 
nitrobenzene by the addition of nitric acid, it was found 
experimentally that, when three equivalents of nitric acid 
were used to one of nitrobenzene at the start, the nitro- 
benzene was half used up at the end of twenty minutes, 
and that, at that time, the ortho-, meta-, and para- forms 
of the dinitrobenzene were present in the proportions of 
6.4, 93.5, and 0.1 parts, respectively. The reaction for the 
formation of each of these three substances is of the second 
order, the rate being proportional both to the concentration 
of the nitrobenzene and to the concentration of the nitric 
acid. Required to find the three velocity constants. 

Let c = no. equivalents nitrobenzene left at time 6 

c' = no. equivalents nitric acid left at time 
x,y,z = no. equivalents ortho-, meta-, para-products 
at time 0. 

The set of three simultaneous differential equations will hold, 

dx 7 , dy j . dz 7 , /r _, 

de = l ' dB = 2 ' d6 = ^ ^ 

Furthermore one equivalent of nitric acid is expended for 
each equivalent of either of the three products formed; 
and the same is true of the nitrobenzene. Hence 

c = c (x + y + z), and c' = c</ (x + y + z) (58) 

The solution may be carried out much as in Ex. (I). By 
dividing one of the differential equations by another and 
integrating it is easy to see that 

x : y : z = fci : fe : fc 3 = 6.4 : 93.5 : 0.1 (59) 

Again, by adding the three differential equations (57) and 
putting x + y + z = u there results 

- (fci + h + W (co - u) (co' - i*) (60) 
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so that clearly the sum fci + kz + kz can be determined by 
an equation of the form (41), because (60) is of the same 
form as (39). The calculation is left to the student. 

Example (K). A certain tank is filled with water. At 
a given instant two orifices in the side of the tank are opened 
to discharge the water. The water at the start is 10 feet 
deep and one orifice is 6 feet below the top while the other 
one is 8 feet below the top. The coefficient of discharge of 
each orifice is known to be 0.61. The tank is a right circular 
cylinder with vertical axis, and is 6 feet in diameter. The 
upper and lower orifices are 2 inches and 4 inches in diameter, 
respectively. How long will be required for the tank to be 
drained down to a depth of 5 feet? 

Solution. Call D the depth of water in the tank at any 
time 6. The head, in feet of water, above the upper orifice 
will be D-4, while the head above the lower orifice will be 
D-2. The discharge from an orifice of this type is known 
from hydraulics to satisfy the law 

u = cV2lfi (61) 

where u is the velocity of discharge through the orifice, (that 
is u is the speed with which a particle of water is moving as 
it leaves. The constant c is the coefficient of discharge, 
which in an ideal fluid free from friction and surface tension 
might be taken as unity in an ideal orifice, but actually is 
found by experiment to be about 0.61 for ordinary small 
orifices with sharp edges; g is the acceleration of gravity, 
32.2, and h is the head above the orifice in feet. The for- 
mula is a statement of the fact that the kinetic energy of the 
water on leaving is equal to what it would have had if it 
had fallen freely through the height h, reduced by departure 
from ideal conditions. 

Call the velocity through the upper orifice HI and that 
through the lower i^. Then 

MI = 0.61 V2 (32.2) (D-4) 
and MS = 0.61 V2 (32.2) (D-2) (62) 
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The volume of water discharged per second through the 
two orifices will be u\Ai and u^A 2j respectively, where AI 
and A 2 are the areas of the orifices in square feet. 

The change in the depth of water in the tank is equal to 
the volume of water leaving the tank divided by the cross- 
sectional area of the tank, therefore 

dD __ UiAi + u> 2 A 2 . . 

de~ 9 TT (M) 

where AI = Tr/144 and A 2 = T/36. Collecting values and 
simplifying 

dD Q.61 VeO 



de 9(144) 
whence the time for D to fall from 10 to 5 feet is determined 

by 

- r . _ dD / _ = 0.00378 C $ dO (65) 
J 10 VD-4 + 4 VD-2 Jo 

The integral on the left is evaluated by first rationalizing 
the denominator, (that is, multiplying above and below by 
A/D-4 4 VD-2), and then separating into the sum of 
two integrals, 



The calculation may now be completed by the aid of the 
Hudson-Lipka table of integrals, or by Simpson's Rule. 

14. Opposing Processes. A second very important 
type of problem arises when two simultaneous processes, 
instead of occurring independently, tend to oppose or 
reverse each other. 

Example (L). A substance C is changed into a second 
substance X by a reaction of the first order. At the same 
time X is changing back into C, also by a reaction of the 
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first order. Observed data: 

When 0=0 5 

C = 10 7 3 

X = 3 7 

Find the velocity constants k L and & 2 of the two reactions. 

Solution. Let c and x denote the equivalents of C and X 
at time 6. An equation like (18) stated with respect to X 
is 

dx = fcic dd & 2 z d0 (67) 

This equation contains three variables, but 

x + c = x + CQ (68) 

because tne total number of equivalents of both substances 
is not altered by either reaction. Eliminating c 

fir 

|g = fc (c + X G ) - (fe + fa) x (69) 

This differential equation has two important properties. 
First, it is evident from the form of the equation that, as x 

dx 
increases, the rate -r- will approach zero. In other words, 

the system will approach a condition of equilibrium as time 
goes on, a well-known fact. 

In the second place, if we adopt the notation (see the 
remark below) 

*" = fcTTV (co + a;o) (70) 

and also 

y = - x (71) 

the differential equation (69) becomes 

) M (72) 



which means, in words, that the rate of the reaction is 
proportional to the departure of the system from a state of 
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equilibrium, and that the process is of the first order with 
reference to the variable y by which this departure is meas- 
ured. 

That a physical or chemical system may be regarded as 
activated by the divergence between the state of the system 
and a state of equilibrium is a principle of wide application. 

To find the velocity constants, first find their ratio, by 
using the observed values in (70). Then find their sum in 
the usual way from (72). The calculation is left to the 
student. 

Remark on the above solution. The quantity x* defined 
by (70) is the particular value of x which would cause 
the right member of (09) to vanish. It is therefore the 
value which x actually approaches as time goes on. It is 
also called the equilibrium value of x. Similarly in many of 
the problems which follow, if a variable approaches a limit 
as time goes on, the subscript will be used to indicate this 
limiting value. 

In the above example, by eliminating the variable c, 

dx 

wo obtained a differential equation (69) which expressed -rr 

au 

as a polynomial of the first degree in x, and it appeared in 
(72) that this polynomial was merely a constant times the 
factor x* x, here called y. 

In general -rr will be found to be equal to some more 
au 

complex expression. In many cases which are important 
in practice, serious algebraic troubles may be avoided by 
using a principle of which the above is a simple illustration: 

cJy 

Theorem. If -,~ is equal to a polynomial in the single 
a(J 

variable x, (together with constants), and if x approaches a 
limiting value x*. as time goes on, then x x is a factor of 
the polynomial. 

This theorem is an application of what is known in al- 
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gebra as the factor theorem, proved in all good text-books of 
algebra. 

In experimental work on reactions or other processes which 
tend to approach a state of equilibrium, it is important, 
when possible, to determine by actual measurement, the 
limiting or equilibrium value of the variable (or variables) 
on which this state depends. If the experiment is stopped 
a long way from this final state, it is always difficult, and 
often quite impossible, to determine accurately velocity 
constants and related numbers. It may even bo impossible 
to decide on the order of the reaction or the nature of tho 
process. The truth of these statements will appear from a, 
study of the problems which follow. 

Example (AT). Opposing reactions of the second order. 

Acetic acid and alcohol react to form acetic ester and 
water, 

IICaHaO* + 0>II 5 OH ^ CWaHaOa + H 2 ^ (73) 

being of the second order both ways. Observed were 
When = 64 oc days 

acid = 1 0.750 0.333 gm. mols. 
alcohol = 1 0.750 0.333 " 
ester = 0.250 0.667 " 
water = 0.250 0.667 " " 

Find how many days arc needed for the alcohol to fall to 
half the original amount. 

Solution. We shall as usual first set up the general 
equation for such a process, and obtain a result in which 
any particular data may be substituted after the integration. 
Adopt the notation 

x, y, u, and v = gm. mols. alcohol, acid, ester, and water 
at time 6, respectively. 

Then the differential equation for the alcohol is 



i (74) 
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To eliminate y, u, and v use the stochiometric relations 

XQ = 3/o y = U UQ = V VQ (75) 

which are equivalent to 

y = X + t/o XQ, U = XQ + UQ X, V = XQ + VQ X (76) 

and by virtue of which the differential equation (74) becomes 
^5 = ax 2 + bx + c (77) 

where 

C = k 2 (Xo + UQ) (XQ + VQ) (78) 

and are constants for any particular run. Now by definition 

Of Xoc 

= ax^+bx* +c . (79) 

because as x approaches the limiting value # the chemical 

dx 

system approaches a state of equilibrium and -rr approaches 

au 

zero. Subtracting (77) from (79) 

- j~ = a (x 2 <x - z 2 ) + b (x - x) (80) 

the constant c cancelling by subtraction. The right side of 
(80) factors easily and we have 

-~= (zee - x) (a*. + ax + b) (81) 

in agreement with the factor theorem. It is best, for reasons 
soon to appear, to factor out the constant a also, giving 

_ ^ = a t x _ x \( x I x I \ (Q<>\ 

1 s\ I* ^t/ CC Jj J I ^ I *' OC I I \^~/ 

(I/U \ U/ 

Much the easiest way, in practice, to handle this differential 
equation is to notice first that, by setting the right side of 
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(74) equal to zero (when of course the variables take their 
equilibrium values), the condition for equilibrium may be 
written in the form 



The importance of this equation resides in the fact that, if 
we know the limiting values of the variables, we thereby 
know the ratio of the velocity constants, and furthermore, 
even if unfortunately we are unable to observe by experi- 
ment any one of these limiting values, we shall still make 
our calculations much simpler by introducing this ratio as 
one of our unknowns. We shall accordingly represent this 
ratio by the letter m, namely 

- (84) 

We now return to the fraction - occurring in the second 

factor on the right of (82). By using the values of a and b 
from (78) and introducing m by (84) we find 

6 1 



a m I 



[(1 2 ra) XQ 7/0 rnuo - mv {) ] (85) 



By inspection of this result it appears that if the equilibrium 
values of the variables be known as well as their initial 

values, then - is known and (82) may be written 

fir 
- 5 = o (*-*)(* + 0) (86) 

where 

P = x+^ (87) 

The integration of (86) may now be carried out and the value 
of a determined in a manner very similar to Ex. (G), for 
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in form (86) resembles the differential equation of a simple 
second order process. Using the data above given for 
alcohol and acid we find ki = 4 & 2 or m = J. The differ- 
ential equation (80) reduces to 

air^r ($-*)(* + !> (88) 

The completion of the work is left as a problem for the 
student. 

It may happen in some special case that ki = k. The 
constant a is then zero, and cannot, be factored out. But 
(81) reduces to the form of a first order process and presents 
no difficulty. 

Example ( N). Reactions of the first order in one direction 
but of the second order in the reverse direction. 

Carbon monoxide decomposes into carbon and carbon 
dioxide in the presence of a large excess of carbon. The 
equation is 

2 CO ^ C + C0 2 (89) 

The rate of decomposition of the monoxide is proportional 
to the square of its concentration, while the reverse re- 
action proceeds in proportion to the first power of the 
concentration of the dioxide. These reaction rates are 
commonly explained, very roughly, by saying that, under 
the kinetic theory of gases, the rate of any reaction depends 
on the probability of a collision between the molecules 
concerned. A more rigorous explanation can be deduced 
from the laws of thermodynamics, by Gibbs' theory of 
chemical potentials. 

The above statements as to the rate of the reactions are 
under the assumption that the reaction takes place at con- 
stant temperature and constant total volume. 

Let x and y be the concentration of CO and C0 2 respec- 
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lively, at time 6, expressed cither in equivalents per unit 
volume, or in partial pressure as preferred. Then 

^ = fci** - ky (90) 

Now the volume, (or pressure), of CO used up is twice that 
of the C0 2 formed from it, hence 

#o - x = 2 (y - 2/ ) (91) 

and eliminating y 

g=/c,^+|x-|(x + 2 ? /o) (92) 

while at equilibrium 

= fc^cc + ^ * - ~ (x + 2 2/ ) (93) 

so that by subtraction 

-g-fe (,-*) (* + *.+-*!.) (94) 

If we set 

m=|, |8-*+f (95) 

we shall have by (90) 

m = ^ (96) 

whence m is known from the equilibrium conditions, on 
which, for most gas reactions, better experimental work has 
been done than on the rate of reaction. The differential 
equation takes the form 

^ = 2fc (* - *)(*+# (97> 

and presents no new difficulty. 
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Example (0). Reactions of the second order one way 
and of the third order in the reverse way. There are a 
number of gas reactions of this type, for instance 

2 CO + O 2 = 2 CO 2 

2 SO 2 + O 2 = 2 SO 3 

2 H 2 + 2 = 2 H 2 

Using the first one for an example, (the theory is the same 
for all), let the concentrations be x, y y and z, for monoxide, 
oxygen, and dioxide, respectively. Then 



(98) 
By setting the right side equal to zero 

-*?? > 

We have also the stochiometric relations 

x x Q = z z = 2 (y 7/0) (100) 

whence, eliminating y and z and introducing x* in the same 
manner as in the preceding example, 

^ = ?(*- -x)(x* + bx + c) (101) 

where 

6 = x X Q + 2 (T/O m), c = bx* + 4 ra (x + o) (102) 

If the initial conditions, as well as the equilibrium conditions, 
are known, b and c are numbers which can be easily com- 
puted. The variables may now be separated, 

(103) 



(xc, x) (x 2 + bx + c) 2 
The left member is the same as 



_ 
x x 2 + bx + c 
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where 

1 *\Ar 

A = B = -fri- T- and C = ^ (105) 

a; 2 * + &Zcc + c a;* 

whence A, B, and C are at once calculated and the in- 
tegration carried out by the aid of the tables. 

(For the algebra, see the " method of partial fractions " 
in any good text book of advanced algebra, or of integral 
calculus, for example, Hall and Knight's " Higher Al- 
gebra " Chap. XXIII. It is oasy to verify that, in general, 
re* is not a root of the quadratic x 2 + bx + c. If such 
were the case, the separation into partial fractions would 
be slightly modified.) 

Any other cases of reversible reactions may be treated in 
a similar manner, provided the initial and equilibrium con- 
ditions are known. 

In many processes other than chemical reactions, a vari- 
able tends to diminish for one cause and to increase for an- 
other cause. The differential equation for this variable is 
often similar in form to that of a reversible reaction. The 
following is a simple illustration. 

Example (P). Suppose the bottom of the tank in ex- 
ample (C) to be covered with a cake of undissolved salt. 
Assume the surface exposed to be substantially constant, 
and the salt to dissolve at a rate proportional to the differ- 
ence between the concentration of the solution and that of 
a saturated solution, (3 Ibs. per gal.). With all other con- 
ditions as in example (C), express the total amount of 
dissolved salt as a function of the time. Given that, if 
the water in the tank were fresh, one pound of salt would 
dissolve per minute. 

Solution. Equation (18) stated for the dissolved salt is 
dc = k (300-c) d8 - 0.03 c dO (106) 

dc 

We have given that if c were zero then -^ would be 1 Ib. per 

du 
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min. Hence fc = ^7^ The differential equation may now 
oUU 

be written 



dO ~ 30 
whence 

c-30 = 20 e~ Q W (108) 

showing that the amount of dissolved salt approached 30 
Ibs. as time goes on (assuming excess of solid salt). 

In the preceding examples the data have been, or have 
been assumed to be, such that the limit approached by the 
dependent variable could be directly found. If not, the 
problem is harder, first by reason of the purely mathematical 
difficulties which arise, second because the data themselves 
may not suffice. The following example will illustrate 
one possible method of attack. 

Example (Q). Diphenyl-ehlor-rnethan reacts in ethyl 
alcohol yielding hydrochloric acid and an ester. The re- 
action is assumed to be of the form 

R01 + R'OH ^ HC1 + ROR' (109) 

It is desired to verify this assumption, and to determine 
the velocity constants. The following data were obtained, 

When 6 = 13 119 142 162 182 212 minutes 

HC1 - 0.00346 0.02680 0.0309 0.0343 0.0375 0.0418 formal 

It is evident that equilibrium was not even approximately 
attained. 

Let x = rnols chlor compound at time 0, 

y = mols HC1 and also mols ester at time 6. 

If the reaction is according to (109) the differential equation 
must be 

= k& - fc 22 / 2 (110) 
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where /bi and & 2 are constants, for the amount of alcohol was 
so large, (about 21 mols), as to remain substantially con- 
stant. In this experiment x = 0.09966. By stoichiomctry 
x + y = XQ whence x may be eliminated giving 

& = fc (* - y) - W (111) 

As before, let y<* be the limiting value of y, so that 

= fci (x, - T/.X) - % 2 (112) 

and by subtraction 



(113) 

a/2 / 

We also have directty from (112) 

(114) 



x<> - y 

and may also write 



^ . 
^ ^ 



so that both &i and j3 can be determined if ?/ is known. 
The differential equation (113) now takes the form 



of which the integral may be written 



The method consists in now assuming a value for the 

unknown 7/<r, and plotting the quantity log - - against 

2/ "~ y 

6 on ordinary cross-section paper, when, if the true value 
of i/oc has been selected, and if the hypothesis as to the 
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nature of the reaction is correct, the result should be a 
series of points as nearly on a straight line as the accuracy 
of the experimental work will justify. Fig. 3 shows the 




Experimental Plot for Example (Q) 
FIG. 3. 

results for trial values of y<* equal to 0.070, 0.075, 0.080, 
and 0.083, indicated by the signs +, X, <, and >, respec- 
tively. The first and the last trials indicate a very slight 
concavity upward and downward, respectively although the 
departure from a straight line is in every case comparable 



SIMULTANEOUS PROCESSES 55 

in order of magnitude with the experimental errors. Values 
of fc 2 estimated for the four trials were 0.0170, 0.0118, 0.0082, 
and 0.0064. Thus the first is more than double the last. 

The only conclusion to be drawn with regard to fc 2 from 
such a set of data, (very good in themselves but not carried 
far enough,) is that, if the reaction is of the assumed form, 
the constant fc 2 probably lies within or very near the interval 
from the largest to the smallest value obtained. The 
nearer the experiment is carried to equilibrium, the more 
accurately fc 2 can be estimated by this method, for the 
differential equation shows that fc 2 is more important the 
larger y becomes. 

With ki it is quite another story, and the above data 
suffice to determine its value very well, because ki is the 
predominating constant at the beginning of the reaction. 
In fact all four trials yield nearly the same value of k\. 
The verification is left to the student. 

The above data were obtained by a student who was 
primarily concerned with fti. In a later experiment with 
the same substances to determine fc 2 , a full series of obser- 
vations, extending to equilibrium, was made by Mr. A. A. 
Morton of the Massachusetts Institute of Technology, 
through whose courtesy we are able to publish the results. 
These will be found in problem 16 below. 

15. Consecutive Processes. A third fundamental type 
of process may be viewed as the result of two simpler proces- 
ses taking place in succession. 

Example (R). Suppose a substance of initial amount c 
changes into a first product, which in turn changes into a 
second product, the rates of the two reactions being in 
general different. (Such reactions are very common, the 
changes, for example, of radioactive substances are of this 
character universally, so far as is known.) Suppose both 
reactions to be of the first order. 
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Let c = amount of parent substance at time 6. 

x = amount of first product FORMED up to time 0. 

y = amount of second product formed up to time 6. 
Then x y = amount of first product PRESENT at time 6. 
The differential equations for x and for y will then be 

=Mc -aO and ^ = fe (x - y) (118) 

because, by hypothesis, x increases at a rate proportional 
to the value of c, that is c x, while y increases at a rate 
proportional to the value of the first product actually present, 
that is x y. The first of these differential equations, 
since it contains only one dependent variable, can be at 
once integrated, giving 

x = Co (l - e-**) (119) 

This value of x is now substituted in the second differential 
equation. The result is 

jj + fey = feco (!-*-*') (120) 

a differential equation " linear " in the dependent variable y. 
The easiest way to solve a linear equation is to write it 
in the form 

^ e +Py = Q (121) 

and then to set 

R = ef (122) 

when the solution will be given by 

y R = CliQdd + const. (123) 

as may be easily verified by taking the value of y from (123), 
putting it back in (121), and showing that the two sides 
become identically equal. 
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In the present example P = fc 2 so that R = #* and (123) 
becomes 

(1 - e-W) d0 + const. 



+ . o2 . c fa-*)* + const. (124) 

A/1 rC2 

To determine the constant of integration, we have y = 
when 6 = 0. Hence 



= C + rj- + const. (125) 





A/2 



so that the constant is equal to , ' by simplifying. 

A/] A/2 

Substituting this value of the constant in (124), and dividing 
both sides by e 

'- w -'- M (126) 



an important equation because it shows the manner in 
which the final product of a pair of first order consecutive 
processes builds up. 

It is sometimes possible to make a separate or independent 
study of the two reactions and thus find the velocity con- 
stants fci and & 2 . This is not always convenient. We 
shall now suppose that the parent substance was originally 
pure, and call its amount unity, that is CQ = 1. Suppose 
also that a complete analysis of a sample was made at a 
known time 6 giving at that time an amount 0.4 of the first 
product actually present, 0.3 of the second product, and 
hence of course 0.3 of the parent substance unchanged. 
Suppose 6 3. That is, we have found by actual analysis 
when 6 = 3 that x y = 0.4 and y = 0.3. We shall show 
how the " transcendental " equation (126) may be solved 
for the velocity constants. 
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In the first place, fa may be regarded as known, for the 
first reaction is a simple first order process, and fa is directly 
found by an equation like (7), Chapter I. 

We then must solve (126) for the unknown quantity fc 2 . 
We arrange the equation with the exponential term e~ k # 
alone on one side, at the same time getting rid of the other 
exponential term e ~ kl by putting in its place its value from 
(119). As a result, (126) takes the much simpler form 



e -w = fa . -Z- + 1 - (127) 

Co/Ci Co 

where by hypothesis all the constants are known except 
fa. We then take a small piece of ordinary plotting paper 
and make a graph of the function e ~ k * regarding & 2 as the 
variable. That is, by assigning a series of values to fe>, 
we find e~^ by aid of a table of exponentials and plot 
e -w as ordinate against & 2 as abscissa. The value of 6 is 
kept constant, arid is the value at the time the analysis was 
made. The result is of course a curve of very familiar 
character. 

We next plot the right side of (127) in a similar manner. 
It is obvious that the result is a straight line, so that it is 
sufficient to find two or three points, calculating the right 
side of (127) for ordinates to match selected values of k% as 
abscissa, drawing a line through these points with a ruler. 
It will in general be found that this straight line cuts the 
curve already drawn for e ~ ** in two distinct points. One 
of these, if the graphs have been carefully made, will have 
for its abscissa a value precisely equal to the value of fa 
already found. The other point will have for its abscissa 
the value of fa required in the problem. 

The reason ki also appears as an intersection is that equa- 
tion (127) actually has two solutions, one of which is given 
by setting fc 2 = fa, a fact easy to verify directly from (127). 
The solution fa = fa is " extraneous " and was introduced 
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when (126) was cleared of fractions to obtain (127), because 
it was necessary to multiply by fci fe. This extra point 
of intersection is a help, not a hindrance, because if we 
know that the line and the exponential curve intersect when 
fc 2 = ki we have a check on the draughtsmanship. 

If in a particular case it should really be true that fc 2 = fci, 
the straight line would be precisely tangent to the exponential 
curve. When the value of fc 2 has been approximately found 
from a first graph, a more accurate value may be found by 
plotting on a larger scale a portion of both the curve and the 
line in the required region. Thus fc 2 may be found as closely 
as needed. 

Applying this method to the data above assumed to be 

known, we have first fci = - 5 In (0.3) = 0.40 and then (127) 

o 

becomes 

e -M3) = _ fc 2 + 0.7 (128) 

By plotting, it is easily found that fc 2 is near the value 0.35. 
As a final determination, a few values may be tabulated from 
a table of natural logarithms, thus 

In 0.35 = - 1.04982, fc 2 = - In 0.35 = 0.34994, 

- fc 2 + 0.7 = 0.35006 

In 0.355 = - 1.03564, fc 2 = - I In 0.355 = 0.34516, 

- fc. 2 + 0.7 = 0.35484 

In the first line, we have taken the value of e ~ 3 * 2 to be 
0.35, giving the right side of (128) as 0.35006, which is 
larger than 0.35. 

In the second line, we have assumed the left side to be 
0.355, giving the right side as 0.35484, which is smaller than 
0.355. Hence fc 2 is between 0.34994 and 0.34516. This 
method may be used to check or to replace the graphic 
method, which is illustrated by Fig. 4. 

If fci is not very different from fc 2 , the straight line will be 
nearly, but not quite, tangent to the curve, and will run 
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nearly parallel to the curve near the desired point. We 
may, however, largely overcome this difficulty by greatly 
increasing the scale of measure on the vertical axis, for the 

A 
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0.33 



0.32 



0.31 



t).30 




Axis of k 9 





"6.34 0.35 0.36 0.37 0.38 0.39 

Plot of the Transcendental Equation in Ex. (R) 

FlG. 4. 



0.40 



curvature of the curve is exaggerated while the line remains 
straight. In such cases, however, the method of numerical 
tabulation is to be preferred. In all cases a rougher graph 
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should first be made to determine the approximate value of 
fe. 

Example (S). Show that in a pair of consecutive first 
order processes the first product rises to a maximum and 
then falls off toward zero. Find the maximum amount of 
this first product, and the time when the maximum is reached. 

Solution. Subtracting (126) from (119) and simplifying 



Differentiating with respect to 9 we have 



If we were to plot the graph of x y against 0, equation 
(130) shows that the slope at the start must in all cases be 
Cofci, because the right side reduces to this value when = 0. 
To find the maximum we set the right side equal to zero, 
giving 

|i= e (*i-*>)' (131) 

which solved for gives 

= 7-^-r- In yi (132) 

// IT if x f 

as the time when the first product attains its maximum 
amount. For example, if it were desired to obtain the 
greatest possible yield of this first product, (132) gives the 
proper time to stop the operation. 

To find the value at the maximum, set, in (129), 

- M -E'- M (133) 

because the right side of (130) is zero at the maximum. 
Then (129) simplifies to 
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-k l e (134) 

at the maximum of x y. The value of 6 from (132) is 
now used, and (134) reduces to 



x - y = co (135) 

as the maximum yield of the first product. 

It is evident that the form of the above investigation 
would require modification in case fci is substantially equal 
to & 2 . This case is left as a problem for the student. 

16. Processes of Mixed Type. The three types which 
have been studied in this chapter, side by side, reversible, 
and consecutive, very commonly occur together. The 
number of possible combinations is very large, but methods 
in general similar to the foregoing may be used to set up 
the differential equations, which, however, are not 
necessarily either separable or linear. 

The following simple but important example will serve 
as an illustration. 

Example (T). Steam was passed over heated carbon, 
yielding " water gas/ 7 which is a mixture of steam, carbon 
dioxide, carbon monoxide, and hydrogen. Reckoned on a 
basis of one unit volume of pure steam at the start, the 
following results were obtained: 

When volume of steam = 1.000 0.906 0.709 0.556 0.376 0.056 
volume of CO 2 = . 000 . 041 . 100 . 123 . 099 . 0242 

Show that these results are in good agreement with the 
hypothesis that the reactions are 

(A) C + 2 H 2 O = CO 2 + 2 H 2 

(B) C + H 2 O = CO + H 2 

(C) C + C0 2 = 2 CO 

together with the hypotheses that all three reactions are of 
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the first order and that none of them are reversible to a 
perceptible extent under the conditions of the experiment. 

Also determine the RATIOS of the velocity constants. 

Solution. Let fc b & 2 , and fc 3 be the velocity constants of 
the respective reactions. The carbon dioxide is formed in 
(A) and used up in (C). Letting x and y denote the volumes 
of steam and of C0 2 at time 6, the differential equation for 
y will be, on the above hypotheses, 

=*i*- ^y (136) 

Again, the steam is used up in both (A) and (B), these 
two being side reactions, while (C) is consecutive. In 
writing (136) the constant fci was chosen to correspond to 
the rate of formation of one volume of C0 2 , corresponding 
to two volumes of steam. The differential equation for 
the steam will thus be 

flr 

j= - 2 fc,z - Jfcax (137) 

Since it was not practicable to observe the values of 0, 
we divide the two differential equations to eliminate dO. 
The result is 

dy _ kix fa?/ 
dx 2 k\x k 2 x 

If we divide numerator and denominator by fci and, for 
simplicity, adopt the notation 

=^ and a = 2+^ (139) 

K,\ KI 

the equation (138) takes the form 

dy __ $y - x 

dx - ~^T 

This equation is linear in y and may be integrated as in 
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the previous example. The values of a and may then be 
found by trial, hence the ratios fc 3 /fci and fe/fci are known. 
The computation is left to the student. 

Remark on the above solution. We might perhaps expect 
reaction (A) to be of -the second order, a hypothesis which 
leads to a differential equation not in agreement with the 
data. We must remember, however, that this reaction 
takes place on the surface of the carbon, so that we should 
not reason by analogy with gas reactions. For a full dis- 
cussion, with many more data, see a paper by Haslam, 
Hitchcock, and Rudow, in " Industrial and Engineering 
Chemistry," Feb., 1923. It is of both theoretical and prac- 
tical interest that data over a wide range of pressure, and 
also of temperature, agree with the equation (140). 

PROBLEMS 

1. A substance c yields two products x arid y. The following ob- 
servations were made: 

When 0=0 5 

c = 10 7 

x = 2 

y = 1 

Find the velocity constants, assuming two side reactions of the first 
order. Docs it matter in this problem whether the substances are 
measured in grams or in gm. mols? 

2. Complete the calculations in Ex. (J). 

3. Complete the calculation in Ex. (K). 

4. Set up and integrate the differentia! equations for a pair of side 
reactions, one of which is of the first order, the other of the second. 
Does Wegscheider's principle hold? 

6. Complete the calculations in Ex. (L). 

6. It was shown by P. Henry that gamma-oxy-butyric acid is con- 
verted into gamma-butyrolactone at the following rate: 

When 5 = 50 100 220 
lactone = 4.96 8.11 11.55 13.28 

The initial amount of acid was 18.23. Show that these data are in 
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agreement with the hypothesis that the reaction is reversible and of 
the first order both ways. 

Hint: plot x* x against 6 on semi-log paper. 

7. Find the velocity constants in Prob. 6. 

8. Complete the calculations in Ex. (M), finding also the values 
of the velocity constants. 

9. How long would it have taken in Ex. (M) for the alcohol to fall 
to half its value if we had started with three equivalents of acid to 
one of alcohol, other conditions remaining such that the velocity con- 
stants are the same as before? Find also the amounts of the four sub- 
stances present when equilibrium is reached. 

10. By experiment it is known that when carbon, carbon monoxide, 
and carbon dioxide are in equilibrium at 900 C. the partial pressure 
of the dioxide is 0.0232 times the square of the partial pressure of the 
monoxide. Assuming that the experiment started with pure CO in 
presence of excess of carbon, and adopting as unit of time the interval 
required for the CO to diminish by one tenth of one per cent, find the 
time required for the CO to fall to the value midway between its initial 
value and its equilibrium value. Assume one unit of pressure at the 
start. 

11. With the data of Ex. 10, assume that the experiment started 
with pure C02 in presence of excess of carbon. Adopt as unit of time 
the interval required for the CO 2 to diminish by one tenth of one per 
cent. Find the time required for the COz to fall to the value midway 
between its initial value and its equilibrium value. Assume one unit 
of pressure at the start. 

12. What is the ratio of the time units employed in Probs. 10 and 
11? Are the velocity constants the same in these two problems? 
Explain. 

13. Find the per cent analysis of the equilibrium gases in Probs. 10 
and 11. Would the number 0.0232 have to be altered if we altered the 
unit of pressure? Do you know the unit of pressure actually em- 
ployed in this experiment? 

14. With all other conditions as in Ex. (P), assume the water in the 
tank to be fresh at the start. Express the total amount of dissolved 
salt as a function of the time. What limit is approached as time goes 
on? Compare your answer with that of Ex. (P). 

16. A square tank is originally filled with water to a depth of 6 ft., 
but an orifice two ins. in diameter is opened in the bottom. To make 
up the loss, water is let into the top of the tank at a constant rate, 
after which it is observed that the water level falls to four ft. and re- 
mains constant. How fast is the water coming in at the top? 
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16. Diphenyl-chlor-methan reacts with excess of alcohol as in Ex. 
(Q). The following data were obtained: 

The chlor compound at the start was 0.08808 formal. 

When o=l 3 5 8 10 16 

HC1 = 0.001520 0.001979 0.002455 0.003098 0.00367 0.00484 

When 0=21 26 31 41 51 61 

HC1 = 0.00586 0.00693 0.00795 0.00994 0.01190 0.01387 

When = 63 82 103 125 151 175 

HC1 = 0.01430 0.01860 0.02135 0.02492 0.02855 0.03246 

When = 1037 <x minutes 

HC1 = 0.07179 0.07405 formal 

The value of the HC1 f or = <* was obtained by allowing the experi- 
ment to proceed until the HC1 remained constant at 0.07405 for two 
days in succession. 

Show that this run is consistent with the hypothesis that the re- 
action is of the first order in one direction, and of the second order in 
the reverse direction, as in Ex. (Q). Find the velocity constants. 

17. From the data of Ex. (Q) find fa. Compare with Prob. 16. 
Should the values of fa from Ex. (Q) and from Prob. 16 be expected to 
agree? How ought they to be compared? If the alcohol had not been 
in excess what changes would you make in the method of calculation? 

18. Study the following runs with the object, first, of finding whether 
the reaction has the form and character assumed in Ex. (Q), and second, 
of finding the velocity constants. 

Run (a). P-methyl diphenyl chlor methan in ethyl alcohol. 
Initial solution of chlor-compound was 0.09690 formal. 

When = 10 14 18 22 30 41 52 68 
HC1 = 0.0348 0.0439 0.0503 0.0547 0.0606 0.0642 0.0659 0.0668 

Run (6). O-chlor diphenyl chlor methan in ethyl alcohol. 
Initial solution of chlor-compound was 0.09730 formal. 

When - 18 49 93 159 209 251 339 

HC1 = 0.00014 0.00022 0.00032 0.00048 0.00061 0.00071 0.00091 

When e = 400 953 minutes 
HC1 = 0.00105 0.00229 formal 

Run (c). Diphenyl chlor methan and isopropyl alcohol. 
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Initial solution of chlor-compound 0.09970 formal 

When = 45 74 106 134 165 225 958 minutes 
HC1 = 0.0012 0.0020 0.0028 0.0036 0.0044 0.0059 0.0212 formal 

Run (d). P-phenyl diphenyl chlor me than in ethyl alcohol. 
Initial solution 0.02463 formal. 

When0 = 11 21 30 43 55 77 94 125 

HC1 = 0.0075 0.012 0.015 0.0179 0.0195 0.0212 0.0217 0.0222 

19. Find the velocity constants in Prob. 1, assuming that the re- 
actions are consecutive instead of side by side. Assume x to be the 
first product and y to be formed from x, both reactions being of the 
first order, and not reversible. 

Caution. Note that x, the amount of the first product PRESENT 
corresponds to the quantity x y of Ex. (R). 

20. Find the velocity constants in Prob. 1, assuming that the re- 
actions are consecutive instead of side by side. Assume y to be the 
first product, and x to be formed from ?/, both reactions being of the 
first order and not reversible. 

Caution. Note that ?/, the amount of the first product PRESENT 
corresponds to the quantity x y of Ex. (R), and that x } the amount 
of the second product present, corresponds to ?/ of example (R). 

21. Using the results of Prob. 1, plot two curves showing the amount 
of each product as ordinate against time as abscissa. 

22. Using the results of Prob. 19, plot two curves showing the 
amount of each product present at time 0. 

23. Using the results of Prob. 20, plot two curves showing the 
amount of each product present at t ; me 0. 

Compare carefully the graphs of Probs. 21, 22, 23. Suppose we are 
able by chemical analysis to determine one product at any time, but 
are unable to determine a second product directly, can we always 
decide whether the reactions are consecutive or side-by-side? If 
consecutive, can we decide which is the first product? 

21. Using the results of Prob. 19, calculate accurately the time when 
x is a maximum, and find the value of this maximum. Check by the 
graph of Prob. 22. 

25. Uranium changes into radium by a first order process. Radium 
also decomposes by a first order process. If we started with a given 
amount of pure radium, this amount would fall to half its value in 
about 1300 years. In ancient rocks the ratio of the amount of uranium 
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X) that of radium is 3,200,000 to one. How long a time would be needed 
'or a given amount of uranium to fall to half its value? 

Some rocks are actually well over 50% uranium. State carefully 
ust what can be deduced as to the age of these rocks, on the basis of 
;he above facts. 

26. Determine the ratios of the velocity constants in Ex. (T). 

27. Steam was passed over heated carbon, yielding water gas, the 
lempcrature being much higher than in Ex. (T). The following data 
vere observed, 

Mien vol. of steam = 1.000 0.788 0.725 0.532 0.340 0.255 0.110 
vol. of CO 2 = 0.000 0.061 0.065 0.126 0.093 0.067 0.021 
>how that these results are in approximate agreement with the as- 
umptions of Ex. (T) and find the ratios of the velocity constants. Is 
lie slight departure from the assumptions of Ex. (T) progressive or 
incidental? If not accidental, on what hypothesis may it be ex- 
)lained? 

28. Mustard gas (M) suffers hydrolysis in presence of excess of 
vuter, yielding an intermediate compound (I) and hydrochloric acid 
vhich is of course dissociated into ions H + and Cl~. The intermediate 
lornpound reacts with more water yielding a final product (G) and more 
lydrochloric acid. The reactions may be written 

(M) + H 2 o = (i) + H+ + cr 
(i) + n 2 o = (G) + K+ + cr 

mt it is conceivable that cither one of these reactions may be re- 
'ersible. 

The following further facts are known: 

1. As time goes on, the mustard gas is completely converted into 
G), the amounts of both (M) and (I) approaching zero as a limit. 

2. The net rate of hydrolysis is retarded by the addition of acids. 

3. The second reaction under all circumstances proceeds very 
apidly in the direction of the formation of (G). 

With dependent variables (M), (I), (G), (H+) and (Cl~), set up a 
>air of differential equations corresponding to the pair of reactions. 
A^rite down also three stoichiometric relations by means of which three 
>f these variables can be eliminated. Which variables are best to 
sliminate? 

The resulting pair of simultaneous differential equations need not be- 
ntegrated. 



CHAPTER V 
EQUATIONS OF FLOW 

17. Continuous and Intermittent Processes. In the pre- 
ceding chapters have occurred a number of problems on 
flow. While typical and important, these have been math- 
ematically of a fairly simple character, because the concen- 
trations of the substances concerned were assumed to be 
constant throughout the tank, gas-flask, or recitation room 
at a given instant of time. The flow was taken to be slow 
enough so that the incoming current could mix thoroughly 
with the material present. The differential equation of 
the process was therefore concerned only with the time- 
change of concentration in the system as a whole. 

There is another type of process which is just the opposite: 
the concentration of any material is constant at a fixed 
point in the system, in the sense that at that point the 
conditions do not change at all with time, but as we travel 
from point to point in the apparatus we find a change of 
condition. A process of this kind is known as " continuous " 
because it can be carried on for an indefinite period without 
interruption. The state of the system is known as " steady," 
implying no time-change at any given point. The differ- 
ential equations corresponding to a continuous process will 
be characterized by the absence of partial derivatives with 
respect to time taken at a fixed point. The independent 
variable will usually be the distance a particle of material 
has traveled since 'entering the apparatus, or some quantity 
of similar nature defining position in the system. 

A typical example of what is meant by a continuous 
process is afforded by Example (V), below, where the wet 
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clay is fed into the upper end of a cylinder at a substantially 
constant rate, while the dry clay falls out at the lower end, 
the water being carried away by a steady current of air 
flowing in the opposite direction. The conditions at any 
one point in the drier are constant, that is, the state is 
steady. The mathematics of a steady state is usually not 
difficult. In the case of the drier just mentioned, we follow 
the history of a definite portion of clay in its journey through 
the drier, using as independent variable either the distance 
it has traveled since entering, or the time it has been on the 
way, as we may prefer. For dependent variables we shall 
have such quantities as the humidity of the clay, its specific 
gravity, the humidity of the air in contact with it, and so on. 
If we remained at a fixed point in the drier, we should ob- 
serve no changes at all or, as above remarked, time deriva- 
tives at that point would be zero, because the state is 
steady. 

In processes of the most general type, one or more of 
the variables will change both in time and in space. A 
process which is not continuous is said to be " intermittent." 
If clay were placed in a loft, and a current of air were driven 
over it until all moisture had been extracted, the clay being 
then replaced by a fresh batch, the process would be inter- 
mittent. It would clearly not be a continuous process with 
respect to the clay, nor yet with respect to the air current, 
which would gather moisture more rapidly in the early 
stages. 

If the state of some variable is " steady," there will be no 
time derivatives of that variable taken at a fixed point, 
or to say the same thing in a more mathematical way, the 
partial derivative of that variable with respect to time will 
be zero. On the other hand, if we travel along with a portion 
of the material, that variable will in general change, that is, 
its total derivative with respect to time will not be zero. 

The occurrence of partial derivatives makes some of these 
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processes more difficult in a mathematical sense. In many 
cases, however, there are simplifying circumstances which 
make it possible to handle the resulting equations. In 
the present chapter we shall first illustrate the general 
method of setting up the differential equations of a process 
involving flow, and shall afterwards give examples of the 
simplifying circumstances in a variety of cases. 

18. Illustrations of the General Method. The logic em- 
ployed is essentially the same as that of the first and sim- 
plest problem on flow which we had, namely Example (C), 
and consists merely in setting down in mathematical lan- 
guage the fact that matter is not created or destroyed. 
Considering a definite portion of space, this consists in say- 
ing that all material which enters the space must either 
stay there or pass out. More mathematically, the rate at 
which a particular substance increases within that space is 
equal to its rate of income minus its rate of output, which is 
equation (18) under Example (C). 

The portion of space usually considered in setting up the 
equation of continuity (18) is that bounded by two sections 
of the apparatus at the two points x and x + dx, where x 
is the distance from the entrance to the drier, tube, or other 
apparatus. 

Example (/) Assume a compressible gas flowing 
through a long straight tube of uniform cross-section, and 
undergoing a change of density during its passage through 
the tube due to a drop in static pressure or other cause. 

Let 

u = velocity of flow at any point x at any time 6. 
p = density of the gas at point x and time 6. 
A = cross-section of the tube, a constant. 

We may suppose the units employed to be feet, pounds, 
and seconds. 

The mass of gas passing the point x per second is the 
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product of the velocity, the density, and the cross-section 
A, that is Apu. 

The mass of gas passing the point x + dx will be A [pu + 
d (PU)]. 

The difference between these two expressions, that is 
A [d (pu)], will be the rate of accumulation of gas in the 
space between the points x and x + dx. Now if we call 
p' the average density of the gas in this space, the actual 
mass of gas in the space will be Ap'dx, because Adx is the 
volume of the space and mass is volume times density. The 
rate of change of this quantity, that is the rate of accumu- 

lation of gas in this space, is - - (Ap'dx) which is the same as 

ou 

Adx ~r since A is constant and x is independent of 6. 
ou 

Equating the two values just found for the rate of ac- 
cumulation, 



- Ad ( P u) = Adx - (141) 

oa 

We now cancel A, divide by dx, and let dx approach zero. 
The average density p' approaches p as a limit, and the 

quotient ^ > as dx approaches zero, approaches the par- 

r\ f \ 

tial derivative - > because the distance dx was an in- 
dx 

crement in x taken independently of the time 6. The 
above equation thus becomes 

U ' 



dx dO 

which is the equation of continuity for this case. 

In obtaining this partial differential equation no assump- 
tion whatever was made with reference to the character of 
the flow of the gas through the tube. The process may be 
either continuous or intermittent. The equation is a 
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necessary condition for the flow of a compressible gas through 
a tube of constant cross-section, under the sole limitation 
that all of the gas which enters the tube must either remain 
in the tube or pass out, in other words, it is not destroyed 
by any chemical process within the tube, or absorbed by 
the walls of the tube, or in any other way got rid of. 

It need hardly be remarked that such an equation is not 
by itself sufficient to enable us to find the dependent vari- 
ables u and p at any point and at any time. To do so, two 
further pieces of information would be required. 

1. Since there are two dependent variables u and p, 
we shall require, besides equation (142), one further relation 
of some kind satisfied by u or p or both. This relation may 
or may not consist of a second differential equation. 

2. In general, in order to solve completely a partial 
differential equation, such as (142), we need to know the 
values of the dependent variables u and p throughout the 
system for at least one value of the time, (as at the start of 
the experiment), and need furthermore to know their values 
at some one point (as at the entrance to the tube) for all 
values of the time. This statement will naturally require 
amplification or modification in special cases, and the in- 
formation here described may frequently be replaced by 
other facts, which, however, would serve a similar purpose. 

As a special case, we may now suppose the process to be 
" continuous," which is the same as assuming all partial 
derivatives with respect to time to be zero. The equation 

(142) becomes 

d (pu) 



dx 



= (143) 



or merely J* - = 0, because we have now only one inde- 
pendent variable. Hence 

pu = constant (144) 

a simple but important relation for a process of this type. 
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In this discussion we have used the words " compressible 
gas," but there is nothing in the reasoning that does not 
apply equally well to any substance whatever, provided 
that in reckoning density we make the calculation with 
respect to that substance alone, neglecting other material, 
such as moisture, which may have been lost or gained in 
the process. The next example will serve as illustration. 

Example (V). A lumpy material is dried in a hollow 
rotary cylindrical drum 5 feet in diameter and 30 feet long. 
The drum is substantially one third full of lumps at all 
points. The axis of the cylinder is nearly but not quite 
horizontal. The material enters the drum at a constant 
rate and travels slowly through it by gravity due to the 
tumbling action obtained by the rotation. The lumps 
entering contain two thirds water, which amounts to two 
pounds water per pound of clay reckoned on a water free 
basis. The clay on leaving the drier contains 0.1 pound 
of water per pound of clay. 

The volume of the clay may be taken to be a linear func- 
tion of the water it contains. The material entering weighs 
31 pounds per cubic foot and leaving 22 pounds per cubic 
foot. 

The drier discharges 550 pounds of material per hour. 
The rate of evaporation may be taken as proportional to 
the water content of the material at any point. 

Let it be required to set up the equation of continuity for 
the clay and to find how long it takes the material to pass 
through the drier. 

Solution. We note first that the process is of the con- 
tinuous type. At the exit, we are told that 550 pounds of 
material pass out per hour, containing 0.1 pounds of moisture 
per pound of clay. Therefore 500 pounds of clay pass out 
per hour. It follows that, because the state is everywhere 
" steady," 500 pounds of clay per hour pass across any and 
every section of the drier. If this were not so there would 
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be an accumulation or piling up of clay contrary to the 
hypothesis that the state is steady. 

We may next let V represent the volume occupied by one 
pound of clay, reckoned on a water-free basis, (that is of 
one pound of actual clay, not counting the water which it 
contains). Let w be the weight of water held by this 
pound of clay. By hypothesis, V is a linear function of w, 
that is 

V = aw + b (145) 

where a and 6 arc constants. 

Also let u = velocity of material at any point x y and A 
cross-section of material, which is one third the cross-section 

of the drier. The density p of the clay itself will be -^ 

The weight of actual clay passing any point per hour will 
be up A, and, as already seen, is equal to 500, hence 

upA = 500 (146) 

which is a case of equation (144), the equation of con- 
tinuity for the clay. To give this the form of a differential 

dx 
equation, set u = -JT where 6 is the time required for a 

CLu 

portion of clay to move from the entrance to the point x. 
Putting also for p and V their values we obtain 

6) (147) 

This differential equation contains three variables, x, d, 
and w. A second differential equation is yielded by the 
fact that we have assumed the rate of evaporation to be 
proportional to the moisture content, or 

-*- < 148 > 
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which implies that the air current is rapid enough so that 
the humidity of the air in the drier may be taken as constant. 
The values of a and b may be found by using the initial and 
final states of the clay, as given. The pair of simultaneous 
differential equations may be integrated without difficulty, 
and the time of passage through the drier computed. This 
is left as a problem for the student. 

Example (W). As a further illustration of the general 
method of setting up an equation of continuity, let there be 
a cube of solid material of uniform character, and of initial 
constant temperature TI throughout. Let one face of the 
cube be placed in contact with a cooler which is kept at the 
constant temperature TV Let the other five faces of the 
cube be perfectly insulated, so that the heat flowing in or 
out through them may be neglected. Set up the equation 
for the flow of heat in the cube. 

Solution. Consider a cross-section of the cube parallel 
to the first face and distant x from it. Call the specific 
heat of the material s and its density and conductivity 
respectively p and k. Let a denote the length of an edge of 
the cube. By Newton's law of cooling, the amount of heat 
Q that crosses an area A perpendicular to the z-axis in unit 
time is given by 

AT 

Q= -kA^ (149) 

where T is the temperature at the point x. The amount 
which crosses an area at the point x + dx will therefore be 



and the difference between these two expressions, 

(151) 
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will be the rate at which heat is accumulating between the 
two sections, by reasoning quite similar to Example (U). 
Now if T r is the average temperature of the portion of the 
cube between these two sections, the rate of accumulation 
of heat in this portion will be 

spAdx^- (152) 

because Adx is the volume, pAdx is the mass, and hence 
spAdx is the heat capacity of this portion. Heat capacity 
multiplied by rate of change of (average) temperature gives 
rate of heat increase. We now equate our two expressions 
for rate of heat accumulation, cancel A, divide by dx, and 
let dx approach zero, precisely as in Example (U). The 
result is 

7 ^ T dT /1KOX 

k -<& = * p Te (153) 

which is the equation of continuity for one-dimensional 
heat flow. 

This partial differential equation contains only the one 
dependent variable T. We know the initial state, and the 
law governing the conditions on the bounding faces of the 
cube. We shall therefore expect to be able to find the 
temperature distribution in the cube at a point x and time 6. 

The easiest method (when it can be used) for solving a 
partial differential equation with two independent variables, 
is to assume a solution of the form yz, where y is a function 
of one independent variable only and z is a function of the 
other independent variable only. That is, we may write 

T = yz (154) ' 

where we suppose T 7 to be a particular solution of (153), 
with y a function of x and z a function of 6. Substituting 
yz for T in (153) gives 



78 DIFFERENTIAL EQUATIONS IN APPLIED CHEMISTRY 

, d 2 y dZ /^rr\ 

kz = spy (155) 



and by separation of variables 

ld*y Idz 



(156) 



where, for compactness, - has been written for the constant 

quantity - Now the left side of this equation is by hy- 
A/ 

pothesis not a function of 0, and the right side is not a func- 
tion of x, and since the two sides are equal they must both 
be free from x and 0, that is, both sides are equal to a con- 
stant. Let this constant be c 2 , giving 

Id 2 y 9 , 1 dz 9 ,*r*7\ 

- -ri = - c 2 and -^ = - c 2 (157) 
ydx 2 zrdd 

a pair of ordinary differential equations. A solution of the 
second is 

z = -<** (158) 

which shows why we took the constant of the form c 2 , 
for T on the whole will diminish as time goes on so that this 
constant must be negative. A solution of the other equation 
is either 

y = sin ex or y = cos ex (159) 

as can be seen by substitution. To decide between these 
two, we may, without loss of generality choose a temperature 
scale (by writing t = T T 2 ) such that the temperature 
of the cooler is taken as zero. Thus t = when x = 0. 
Then the form sin ex is alone admissible. Collecting re- 
sults, we take as a particular solution of (153) 

t = smcx.e-<** (160) 
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which, it is easy to verify directly, is actually a solution, 
independent of the value of c. It is also clear that this 
function is still a solution if we multiply it by any constant. 
We may next determine the constant c to agree with the 
condition that at the face of the cube opposite the cooler, 
that is when x = a, there is no flow of heat. By Newton's 
law (149) this is the same as saying that the derivative 

vanishes when x = a. From (160) 

^ = ccoscx.e- (161) 

ox 

and since the exponential function cannot be zero we must 
have 

cos ca = (162) 

giving the general solution 

ca = (2 n - 1) I (163) 

where n is an integer. Thus c may have any of an infinite 
number of values. This is a very fortunate circumstance, 
for we have still to satisfy the condition that the temperature 
is constant all over the cube when = 0. Up to this point 
the work of solution of the partial differential equation has 
consisted of a number of very elementary steps of the sim- 
plest character from the point of view of computation. 
There has been nothing in the calculus employed which 
goes beyond the ordinary formulas of the freshman year. 
To satisfy the condition that the temperature is constant 
at the start, while not difficult, is a little less elementary. 
Note first that if n is an integer then 2 n - 1 is an odd in- 
teger. Thus if n = 1, c = - ; if n = 2, c = ^ ; if n = 3, 

& (I 41* 

c = - , and so on. We now assume 
2a 
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t = AI sin Cix . e ~~ C?T * + A 2 sin c^x . e ~ c ^ 

+ 4 3 sin c& . e - <** + - (164) 

an infinite series where the successive values of c are those 
determined by (163). 

Each term of the series, from its method of derivation, 
satisfies the original partial differential equation (153); 
each term vanishes when x = 0; and each term gives a 
derivative which vanishes when x = a. Hence the same 
properties hold for the whole series. When 6 = the 
exponential factors all reduce to unity, and the scries be- 
comes 

Jo = AI sin CiX + A 2 sin c^x + A s sin c& + (165) 

The final step in the solution is to determine the constants 
AI, A< 2 , As, in such a way that to shall equal the con- 
stant T\ T%, so that the initial temperature of the cube 
shall be T\. This is one of the simplest examples of " Fou- 
rier's Series." It is known that 



when < x < w. (Sec References below.) 

whence comparing coefficients between (165) and (166) we 
have 

A l = - (7\ - T 2 ), A 2 = 4- 1 > ^3 = 4- 1 > etc - ( 167 > 

7T o O 

By putting these values in (165) we have as the complete 
solution of the problem 

~" ~i;r- 'e~w + 



\ 

etc.) (168) 

/ 
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by which the temperature of the cube at any point is de- 
termined for any value of 0. 

Remarks on the above solution. This example is one of 
the simplest that can be found on the solution of a partial 
differential equation when the state is not steady. It il- 
lustrates very well the great difference in character between 
ordinary and partial differential equations. The method 
here employed, of building up the complete solution as the 
sum of a series of particular solutions, each of which is a 
product of functions of a single variable, is very important. 
It corresponds to the conception of a musical note as the 
composite of a series of simply related notes known as 
"harmonics"; the resemblance is much more than a 
fanciful one, for vibrating mechanisms can be constructed 
which serve precisely as mechanical models for equations 
very similar to (153). 

For the proof of the expansion (160) see Byerly's " Four- 
ier's Series and Spherical Harmonics," Art. 20. 

For further reading on partial differential equations 
which can be solved by similar methods see, besides 
Byerly, Woods and Bailey's " Analytic Geometry and 
Calculus," Chap. XVII; and Whittaker arid Watson's 
" Modern Analysis," 3rd Kd., Chap. IX. The most recent 
view-point on such problems will be found in " Intcgral- 
Gleichungen," by A. Kneser. 

19. The Problem of Intermittent Extraction. One of 
the most important types of processes in chemical engi- 
neering is the case where one reacting material flows at a 
uniform rate through another reacting material which is 
stationary with respect to the moving material. This is 
most commonly found in intermittent extraction, where a 
solvent is flowing through an extraction chamber containing 
a solid material part of which is soluble in the moving liquid. 
Examples of this may be found in the extraction of tannin 
from bark by means of water, the removal of color from sugar 
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syrup by means of bone char, and the removal of poisonous 
gases from air by gas mask cannisters. 

Consider an extraction chamber of uniform cross-section 
through which liquid is flowing at the uniform rate of u 
pounds per hour. The liquid is removing an extractable 
material from the solid with which the chamber is filled. 
Let 

y = concentration of the extractable material in the liquid, 
z = concentration of the extractable material in the solid. 
a = volume of liquid per foot length of tube. 
ft = volume of solid per foot length of tube. 

Let the symbol d refer to increment of distance x along the 
tube, while A refers to time increments. Consider as usual 
a section of tube between the points x and x + dx. 
Let 

y = average value of y at the point x during time A0. 

Then 

v/wA0 = input of extractable material across the section 

at x during time A0. 
and 
y + dy = average value of y at point x + dx during timeA0 

hence 

(y + dy) uA0 == output of extractable material across the 
section at x + dx during time A0. 

By subtraction 

input minus output = u dy A0 (169) 

which must be equal to the increase in extractable material 
in the space between x and x + dx during time A0. 

Next let y' = average value of y over the space interval 
dx at time 0. Then ay'dx = amount of extractable ma- 
terial contained in the liquid in this space at time 0. 



EQUATIONS OF FLOW 83 

Let z f = average value of z over interval dx at time 0. 
Then ftz'dx = amount of extractable material contained in 
the solid in this space at time 6. Hence 

(ay* + 0z') dx = total extractable material in this space 
at time 6. At the time + A0, the quantity y f has become 
y r + A?/' and z 1 has become z 1 + Az'. Therefore 

[a (y f + Ay') + ft (z f + Az')l dx = total at time 6 + A0. 
By subtraction 

[a.ky r + /3Az'] d# = increase during time A0 (170) 
and by equating values from (169) and (170) 

u dy A0 + ot A?/ dx + ft Az' dx = (171) 

We now divide by dx A0 and let these increments approach 
zero. The average values approach as their limits the 
values at (x, 0). The result is 



This is the equation of continuity for the type of extraction 
process in question. 

PROBLEMS 

1. A current of air, containing 0.0004% C0 2 , is passed through a 
long straight tube of uniform cross-section packed with small lumps of 
freshly prepared CaO. If after passing through 2 ft. of tube the C0 2 
content has fallen to 0.0002%, what will it be after four feet? 

2. With the apparatus of Prob. 1, assume the entering current to 
consist of 50% helium and 50% CO2 by volume. Again suppose the 
per cent CO 2 to fall to half its value after 2 ft. What will it be after 
4ft.? 

Hint. It may still be assumed that the CaO is in excess, but the 
volume changes due to absorption of CO^ are too great to be neglected. 
Assume the rate of absorption of CO 2 at any point to be proportional 
to the concentration of CO 2 at that point. Set up equations of con- 
tinuity both for C0 2 and for He. 
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3. Complete the calculation in Ex. (V). 

4. In the drier of Ex. (V), suppose air enters substantially moisture 
free and leaves carrying 0.1 Ib. water vapor in each Ib. of air (on a dry 
basis). The stock is kept at a temperature such that the water in it 
is in equilibrium with air holding 0.2 Ib. water vapor per Ib. of air. 
The rate of evaporation at each point is proportional to the difference 
between this last figure and the moisture content of the air at that 
point, and also proportional to the water content of the stock. Find 
how long it takes the stock to pass through the drier, all other data 
being as in Ex. (V). 

Hint. St'l up the equation of continuity for the flow of moisture. 

6. Repeat Prob. 4, assuming that air enters, not moisture free, but 
carrying 0.05 Ibs. water vapor per Ib. of air on a dry basis, all other 
data remaining the same. 

6. With the notation of Ex. (W), prove that if the temperature is 
symmetrical about the axis of a cylinder, and does not vary in the 
direction of this axis, 

#7 T 1 oT ^spdT 
(V r r dr k dO 

where r is the distance from the axis. 

7. With the notation of Ex. (W), prove that if the temperature is 
symmetrical about the origin, 

^1 4. ? 2Z! = ^ L 
dr 2 r dr k do 

where r is the distance from the origin. 

8. With the notation of Ex. (W), prove that in general 



= 

dj* dtp & 2 k do 

where #, y, and z are rectangular coordinates. 

9. With the apparatus and data of Prob. 2, set up the differential 
equations of the process, assuming that the CaO is no longer fresh 
enough so that the conditions may be taken as steady. Assume the 
rate of absorbtion of CO2 at any point to be proportional to the product 
of the concentration of CC>2 and the strength of the lime at that point. 
Can you integrate these equations? In these assumptions, what fact 
has been neglected? 

10. In the removal of color from sugar syrup by passing the syrup 
through bone charcoal, the color is adsorbed on the surface of the 
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char under such conditions that, at least during the earl}- part of a run, 
we may assume the equation -.- = ky, with notation of equation (172). 

Furthermore the amount of liquid is so small in proportion to the 
amount of char that the term in a. may be neglected. Integrate equa- 
tion (172) under these conditions and state in words the meaning of 
your result. 

11. With the apparatus and data of Prob. 1, set up the differential 
equations of the process, assuming, as in that problem, that volume 
changes are so slight as to be negligible, (because the concentration 
of COa is small), but taking the CaO as no longer fresh. Assume the 
rate of absorption at any point to be proportional to the difference 
between the concentration of CO 2 in the air and the concent ration of 
OO 2 on the lime. (This is essentially the same as the problem of a 
gas-mask,) Discuss the validity of these assumptions. 

12. With the cube of Ex. (W), suppose a pair of opposite faces to 
be kept for a time at and 100 respectively, the four remaining, faces 
being perfectly insulated. When a steady state has resulted, the hot 
face is suddenly placed in contact with a cooler at 0, both the opposite, 
conducting, faces being maintained at thereafter. Find the tem- 
perature at each point of the cube as a function of the time, given the 
expansion 

/sin x m'n 2 x . sin 3 x 



when < x < TT. 

13. Repeat Prob. 12, if the hot face is left at 100 and the cool face 
is suddenly raised to 100 and kept there. 

14. Suppose the cube to be initially at 100 throughout. A pair of 
opposite faces are suddenly reduced to zero and kept there, the other 
faces being perfectly insulated. Find the temperature at any point as 
a function of the time. 

16. Two iron slabs each 20 cm. thick, one of which is at the tem- 
perature and the other at 100 throughout are placed together face 
to face, and their outer faces are then kept at the temperature 0. 
Find the temperature of a point in their common face and of points 
10 cm. from the common face fifteen minutes after the slabs have been 
put together. Given r = 0.185 in C. G. S. units, and that the series 

? F'^IlL? j_ 2 sin 2 x sin 3 x .sin 5 2 sin 6 x ~| 

is equal to 1 from x to x = ~ an( ^ ls equal to zero from x = = to 

X = rr. 



CHAPTER VI 

THE GRAPHICAL EVALUATION OF INTEGRAL 
EXPRESSIONS 

20. Reasons for Graphic Evaluation. Differential equa- 
tions for physical or chemical processes are often of such a 
nature that integration by elementary algebraic methods is 
impossible. It may also happen that the relation between 
the variables is not known theoretically but only from an em- 
pirical curve derived from experimental data. It is possible 
to evaluate the integrals of such equations graphically. 
It is also frequently possible to save time by evaluating 
graphically integral expressions which can be solved by 
purely algebraic means but which are very time consuming 
when solved in that way. 

Graphical evaluation depends upon the well known fact 
that the integral of the expression ydx is the area under 
the curve obtained by plotting y against x. 

Example (X). In the simple distillation of a mixture 
of two volatile liquids, if the composition of the vapor leaving 
the liquid has been determined experimentally, it is possible 
to derive an equation from which can be calculated the 
composition of the distillate and of the residue after a given 
fraction of the original liquid has been distilled off. 

Lord Rayleigh found the following relations between the 
composition of the liquid and of the vapor for mixtures of 
acetic acid and water, where " composition " means weight 
of acetic acid per unit weight of mixture. 

Suppose that a mixture of acetic acid and water weighing 
W pounds at time 6 and containing a pounds of acetic acid 
and W a pounds of water, be subjected to a simple dis- 

86 
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tillation such that the vapor leaving the still is always in 
equilibrium with the liquid in the still. (Under these 
conditions, the composition of the vapor will have with the 
composition of the liquid the relations of Fig. 5.) The 

composition of the liquid at any time will be ^ with respect 

to the acetic acid, and the composition of the vapor leaving 
the liquid at that time will be found from a curve as 
shown in Fig. 5, constructed from the following data, 
calling the composition of the liquid and the vapor x and 
y respectively. 

Composition of liquid 



1.0 



1.0 



FIG. 5. 



Composition of vapor 
x y 

O.OG77 0.0510 

0.1458 0.1136 

0.2682 0.2035 

0.3746 0.2810 

0.4998 0.3849 

0.6156 0.4907 

0.7227 0.6045 

0.8166 0.7306 

0.9070 0.8622 



Allow dW parts of the liquid to evaporate. This portion 
of vapor will contain da parts of acetic acid. There will 
remain in the liquid W dW parts of a mixture containing 
a da parts of acetic acid, the composition of the distillate 
da 



Since a = xW and y is the composition of the vapor, 

da 
dW 



_ da __ d (xW) _ xdW + Wdx _ Wdx , . 

y "" ^TT7 ~~ x7TT7 "" ,7TT7 "" X ' JW U W 



dW 



whence 



dW 
W 



dx 
y - x 



(174) 
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and by integrating between any two selected values of x, 
which we may call Xi and # 2 , 



(175) 



y - x 



Even though we do not know the equation of the curve in 
Fig. 5 showing the relation between y and x, the evaluation 




i.o 



Fia. 6. 



Fia. 7. 



of the integral on the right side of the equation is practicable 
by the graphical method. Construct a curve with - as 

y ~~ x 

ordinate and with x as abscissa, as in Fig. 6. Measure the 
area under it between the proper limits as shown and place 

it equal to In ^ 

The integral curve of the right hand side of the equation 
may be constructed in the same way, by keeping x\ fixed 
and selecting successively greater values of x 2 and deter- 
mining the area under the curve for each value of x z as above, 
giving a curve of the form shown in Fig. 7, where the ordi- 
nates are the values obtained for the right side of equation 
(175) as described, and are therefore equal to the values of 

In j~ desired. It will be noted that the integral curve is 



asymptotic to x = 1.0 and that if we take a new value for 
Xi we shall get a new curve of similar form. 
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Example (Y). In some cases the plot of y versus x must 
be made with two different scales for x in order to get the 
necessary precision. 

Suppose that the data be plotted as shown in Fig. 8, it 
being necessary to measure the area under the curve between 
the limits of x = 1 and x = 100. In order to get the 
necessary precision in y between the limits of x = I and 
x = 1.5, it is necessary to use the x scale of the size shown 

FIG. 8. 




8 9^100 



1.5 



100 



FIG. 9. 



in Fig. 8. But when this size is used the paper is not long 
enough to include x = 100. Therefore make a new plot of 
the data from x = 1.5 to x = 100 with the same ordinates 
as in Fig. 8 but with the abscissae 1/20 of the ones before, 
giving the curve shown in Fig. 9. 

The area under the curve in Fig. 8 between the limits of 
x = 1 and x = 1.5 may now be measured and found to be 
equal to AI. The area under the curve in Fig. 9 between the 
limits of x = 1.5 and x = 100 may now be measured and 
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FIG. 10. 



found to be equal to A*. Since the abscissae in the second 
figure are 1/20 of those in the first and since the ordinates are 
the same, one square inch in the second figure will equal 
20 square inches in the first. Therefore basing the evalu- 
ation on the first figure, the evaluation of the expression 
desired will be 

Xx = 100 
ydx 
i 

Example (Z). The evaluation 
of definite integrals by the 
graphical method when the units 
of the graph are more compli- 
cated requires considerable care. 
For instance, suppose that it 
was desired to evaluate the in- 

log y dx. This may 



logy 




r* 

tcgral I 
Jx\ 



be done by plotting log y versus 
x and measuring the area" under 
the curve as before. On the 
other hand, the same curve is 
obtained by plotting y versus x 
on semi logarithmic paper. These 
curves are shown in Figs. 10 and 11. The areas will be 
identical only when the distance between log y = and log 
y = 1 on the ordinates of Fig. 10 is the same as the dis- 
tance between 1 and 10 on the ordinates of the semilogar- 
ithmic paper, Fig. 11. 

Example (A A), The problem involved in determining 
the area under a curve between limits of x when the curve 
becomes infinite with respect to ?/, is one frequently met 
especially in connection with thermodynamic equations in 
physical chemistry. 
For instance it may be shown that the minimum work 
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required to completely vaporize a volatile substance from a 
solution containing n mols of that substance at the partial 
pressure p is equal to 



R T 



In dn 
P 




On 2 



where P is the vapor pressure of the pure substance at the 
temperature T, a constant, and R is the gas constant. 
This integral may be evaluated graphically by plotting 

P 

In as ordinate and n as abscissa, measuring the area under 

the curve between the limits of 
n = n and n = 0. The curve has 
the shape shown in Fig. 12. It 
will be noted, however, that p 

approaches zero as n approaches 

p 
zero, whence the value of In 

P 
becomes infinite. The problem 

resolves itself into the determi- 
nation of the area under the 
curve where the curve recedes to infinity at the point n = 0. 

There arc several ways of doing this. One way is to 
measure the area under the curve between the points n\ 
and n% by any suitable means such as by Simpson's rule or 
by means of a planimetcr, selecting the point r? 2 such that the 
laws of dilute solution apply approximately between n 2 
and 0. The remaining area from n^ to may be estimated 
as follows, the area thus obtained being added to the other 
area to give the total value required, 

Assume that between the limits of n z and the value of p 
is a linear function of n in agreement with Henry's law, or 
p = an. Substituting this in the original equation gives 



n 
Fia. 12. 



RT 



fin* 

Jnt an 



dn 



(176) 
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which gives by direct integration 

n,RT(\n + n^i (177) 

\ an 2 / 

which is finite. 

Example (BB). Rittinger's law states that the energy 
needed for crushing a given weight of material is proportional 
to the surface developed. 

In commercial practice when a material is crushed, the 
resulting particles are of varying size, the percentage of 
each size being determined by means of a screen analysis. 
The result of such an analysis is shown in Fig. 13. 



Weight left on 
each screen 




Average diameter of particles on each screen = D 
FIG. 13. 

The amount of energy needed to crush a ton of coal from 
one screen analysis to another may be estimated from Rit- 
tinger's law in the following manner. 

Let w be the weight of coal on any one screen after such an 

analysis, where D is the average diameter of the particles on 

that screen, and let there be n particles of that size. Therefore 

w = an/) 3 (178) 

where a is a proportionality constant. In the same way the 
surface of the coal particles on that screen will be 

8 = 6rcD 2 . (179) 

Since a and 6 are constants, the relation of the surface of the 
particles to the volume will be 

8 = kx/D (180) 

where x is the volume of the particles on the screen, and k is 
another proportionality constant. 
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Assume now that the number of screens is large so that 
only a small number of particles are caught on any one screen. 
The last equation may then be written in differential form 

kdx 



ds = 



D 



The total surface of all the particles on all the screens will 
then be the summation of the surfaces of those on each of 
the individual screens, or 

(182) 



The right hand side of this equation may be evaluated 
graphically, if it is remembered that the volume of the coal 



100 




average diameter of particles on screen 

FIG. 14. 



-1/D 



is proportional to its weight, by taking as a basis 100 parts 
of original coal, and plotting as ordinates the percentage 
which has passed through a given screen diameter D against 
the reciprocal of D as abscissae, as shown in Fig. 14. 




The area of the space under the curve so obtained is there- 
fore proportional to S. If the areas be obtained for three 
screen analyses before and after crushing twice, the ratio of 
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the first two areas will be to the ratio of the second and third 
areas, as the power consumption in the first crushing is to 
that in the second crushing. Such a plot is shown in Fig. 15. 
Example (CC). The vapor pressure of a volatile liquid 
as a function of the temperature is expressed by the Clapeyron 
equation 

d P L 



. 
dT (V - iO T 

which may be simplified by assuming the Perfect Gas law 
to hold and by neglecting the volume of the liquid as com- 
pared with the volume of the vapor, giving the approximate 
form 



_ 
dT RT 2 

where L is the latent heat of vaporization of one mol of the 
liquid and R is the gas constant. 

It is customary to integrate this equation by assuming 
that L is substantially a constant, which is true over short 
temperature ranges. But over longer ones this is not true, 
and great errors may be introduced by such an assumption. 
For instance, the value of L for methyl alcohol varies from 
9.2 at C., to 1.4 at 238.5 C. It is possible to obtain a 
much more precise value for the vapor pressure by allowing 
L to remain a variable and evaluating the integral expression 
graphically. 

The equation may be written 



To evaluate the right side, plot against T 7 , where T 



is the temperature in absolute degrees, R is 1.99 when L is 
expressed in heat units, and the value of L is known em- 
pirically. The area under the curve between the proper 
values of T is the evaluation of the integral expression 
desired. 
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PROBLEMS 

1. Using the data of Ex. (X) plot on graph paper two integral 
curves as illustrated by Fig. 7, using for one of them x\ 0.15 and for 
the other Xi = 0.70. 

fV> e - 0.05 x fa 

2. Evaluate I , ft ni To how many significant figures is 

JQ X -\~ U.U1 

your answer correct? 

Suggestion: Here is a simple case where the "indefinite integral" 
cannot be expressed in finite terms by means of the elementary functions. 
A graphic process similar to Ex. (Y) is easy to use. It is also easy to 
expand the integrand in a series. These two methods may well be 
used to check each other. 

3. Find the value of f hi (1 + x) dx. 



(a) by plotting In (1 + x) on ordinary graph paper, 

(6) by plotting 1 + x on semi-log paper, 

(c) by the usual method of finding the indefinite integral. 

To how many significant figures do your answers check? 

4. Derive (177) from (176). What is the physical meaning of the 
constant a? What data are needed in order to determine the numerical 
value of a? 

6. A certain material was passed through a tube mill at a fixed rate, 
the power consumption being found to be 1.4 K.W. hours per ton of 
material ground. The mill was then doubled in length and similar 
material was fed at the same rate as before. Estimate the power 
consumption required theoretically in the second case from the fol- 
lowing screen analyses. (See next page.) 

6. Discuss the logic of Ex. (BB) by which, in obtaining equation 
(181), a finite quantity s was changed into a differential ds. What 
other steps of similar character must be taken in obtaining (181)? 

7. In Ex. (BB) what would be the value of the ratio b/a if all the 
particles were perfect spheres? What if they were all perfect cubes? 
What if they were all small flat cylinders with diameters ten time their 
thickness? What are we assuming when we take a and 6 to be con- 
stants? 

8. The latent heat of vaporization of methyl alcohol was determined 
by experiment as follows: 

*C. 50 60 70 100 200 230 

cal/gm. mol 9263 8781 8630 8473 7880 4863 2706 

Prom these data plot a smooth curve of L against t . By the method of 
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Ex. (CC) draw a curve from which the rise in the vapor pressure of 
methyl alcohol with temperature can be read off. 

Data for Problem 5 



Mesh 


Opening 
in inches 


Feed 
% left on 
screen 


First 
product 
% left on 
screen 


Second 
product 
% left on 
screen 


2 


.4 









3 


.263 


6 







4 


.185 


10 


2 




6 


.131 


18 


6 




8 


.003 


18 


11 





10 


.065 


16 


17 


7 


14 


.046 


12 


14 


9 


20 


.0328 


9 


14 


12 


28 


.0232 


8 


14 


20 


35 


.0164 


3 


8 


13 


48 


.0116 




8 


13 


65 


.0082 




4 


9 


100 


.0058 




2 


6 


150 


.0041 






5 


200 


.0029 






4 


Through 200 








2 



MISCELLANEOUS PROBLEMS 

The following problems depend upon the principles and methods of 
reasoning illustrated in the foregoing chapters, rather than upon the 
precise formulas which have been obtained in special examples. 

1. In a tank are 100 gals, of brine containing 50 Ibs. of dissolved 
salt. Water runs into the tank at the rate of 3 gals, per min., and the 
mixture runs out at the rate of 2 gals, per min., the concentration being 
kept uniform by stirring. How much salt is in the tank at the end of 
one hour? 

2. In Prob. 1, suppose the outflow passes through a second 100 gal. 
tank initially filled with pure water. How much salt will this tank 
contain at the end of 1 hour? 

3. Suppose that bacteria grow at a rate proportional to the number 
present, but that they produce toxins which destroy them at a rate pro- 
portional to the number of bacteria and to the amount of toxin. Sup- 
pose further that the rate of production of toxin is proportional to the 
number of living bacteria. Show that the number of bacteria increases 
to a maximum and then falls off toward zero, but that under these 
conditions the colony will never become extinct. 

4. Two vertical cylindrical tanks each 10 ft. deep and 4 ft. in diam- 
eter are connected by a short 2-inch pipe at the bottom. If one of the 
tanks is full and the other empty, find the time required to reach the 
same level in both. Assume that the velocity through the pipe is the 
same as that through an orifice under the same effective pressure. 

5. Assuming that the density of sea water under a pressure of p Ibs. 
per sq. in. is 

1 + 0.000003 p 

times its density at the surface, show that the surface of an ocean 5 
miles deep is about 450 ft. lower than it would be if water were incom- 



6. If half the water runs out of a conical funnel in 2 min., find the 
time required to empty. Assume the rate of flow to be proportional 
to the square root of the pressure. 

7. A brick wall is 30 cm. thick. If the inner surface is at 20 C., 
and the outer at C., compute the heat loss per day through 100 
sq. m. of this wall. Take the conductivity as 0.0020 in C. G. S. units. 

97 
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How much coal must be burned to compensate this loss if the heat of 
combustion is 7000 cal/gram and the efficiency of the furnace 60%? 

8. Two steam pipes 20 cm. in diameter protected with 10 cm. thick 
covering of concrete (conductivity 0.0022) and magnesia (conduc- 
tivity 0.00077) are run underneath the soil. If the outer surfaces are 
at 30 C. and the pipes themselves are at 160 C., compute the heat 
losses per sec. per cm. length in the two cases. 

9. A hollow lead sphere (conductivity 0.0827) whose inner and 
outer diameters are 1 cm. and 10 cm. respectively is heated electrically 
by a 10 ohm coil placed within the cavity. What current will keep 
the surfaces at a steady difference of 5 C.? 

10. A wire whose resistance per centimeter length is 0.1 ohm is 
embedded along the axis of a cylindrical cement tube of radii 0.5 cm. 
and 1.0 cm. An electric current of 5 amp. is found to keep a steady 
difference of 125 C. between the inner and outer surfaces. What is the 
conductivity of the cement? 

11. A radioactive substance Ri decays forming R 2 , the period being 
3 weeks. Rz in turn decays to R 3 , the period being three days. (The 
period of a radioactive change is the time needed for the amount of 
decaying substance, if initially pure, to fall to half its initial value.) 
Plot a curve showing the rate at which R* builds up if Ri is initially pure. 
Plot also a curve showing the life of #2. 

12. The natural growth of population would cause the total number 
of inhabitants of the country to double in k years if there were no 
immigration, which, however, brings a constant influx at the rate of 
6 persons per year. Assuming the natural increase among the im- 
migrant population to be the same as among the native element, in 
how many years will the population double? 

13. Repeat Prob. 12, assuming the natural increase among the 
immigrant population to be m times as rapid as among the native el- 
ement. 

14. With the notation and assumptions of Prob. 13, what is the 
maximum value of 6 if the foreign population is never to exceed the 
native? 

15. Repeat Prob. 14 assuming that any person is to be counted as a 
native whose ancestors immigrated p or more years before his birth. 

16. In drying a sheet of heel board, water was lost thus: 

Whenfl = 66 175 251 352 449 491 573 

water = 503.3 384.4 202.8 102.7 50.2 29.6 24.3 16.7 

B = 616 750 minutes 

water = 15.4 15.2 

What was the probable order of the process? 
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17. Hydrogen dioxide H 2 O 2 was observed to decompose in water 
thus: 

Whenfl = 10 20 hours 
H 2 O 2 22.8 13.8 8.25 units 

where the unit of H 2 O 2 is taken as the number of c.c. KMnO 4 solution 
needed to titrate a definite volume of solution. Find the velocity 
constant on the assumption that the reaction is of the first order. 

18. With the data of Prob. 17, find the probable amount of H 2 O 2 
remaining after 40 hours. 

19. In a first order process the following observations were made 
on the product x of the reaction: 

Whenfl = 10 20 30 40 50 60 <x 

x = 60.8 97.7 119.9 133.4 144.4 146.1 153.8 

Determine the value of the velocity constant by the method of averages. 

20. It is empirically true that if the daily production of an oil well 
is plotted on log-log paper against the time since the well started to 
flow, the result will be a straight line. Show that the total production 

up to time is given by the formula P = -O n where a and n are con- 

stants. Show how these constants may be determined from the daily 
production plot. 

21. It does not pay to operate an oil well after the daily flow has 
dropped to an amount 6. With the notation of Prob. 20, show that the 

total yield up to that time is - f - J n- i . Find also the time when 
this yield is attained. 

22. A substance x is formed in agreement with a differential equation 
of the form 



Show how the process may be represented by a straight line on log- 
log paper. 

23. Suppose the amount of the yield x in Prob. 22 to be known by 
means of an empirical curve for all values of 0. (This curve is on ordi- 
nary plotting paper.) Show how to determine the constants k, a, and n. 

24. A substance x is formed in agreement with a differential equation 
of the form 

fif 

= fc (a - Z) (6 - *) 
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and the yield x is known for all values of 0. Show how to determine 
the constants k, a, and b. 

25. Rubber is vulcanized by mixing it with sulphur and other re- 
agents and subjecting the mixture to heat for a period of time. The 
combination of the sulphur with the rubber takes place at a rate which, 
for short periods, is proportional to the amount of free sulphur present 
and also to the difference between the amount of sulphur present as 
combined sulphur and the amount of combined sulphur which would 
be present if all the rubber had combined with sulphur. (100 parts of 
rubber combine with 47 parts of sulphur.) For longer periods it is 
found that the rate of combination is further proportional to a function 
of the time determined empirically to be of the form o^ n + 1 + 0, where 
a and n are constants. Test the following runs with reference to their 
agreement with the above statements: 

Run (a) 
When = 6 8 10 12 14 16 18 20 22 24 min. 

S = 3.70 3.18 3.22 3.17 3.03 2.91 2.84 2.82 2.68 2.53 
where S = parts free sulphur per 100 parts rubber. 

Run (b) 

When = 8 12 16 20 24 28 min. 
S = 1.26 1.00 0.80 0.70 0.50 0.37 

Also estimate the free sulphur at the start. 

26. The rate of solution of a salt in water is proportional to the 
difference in concentration between the saturated film of solution on 
the surface of the crystal and the concentration in the solution. (Where 
the amount of crystals is large, the decrease in their surface may be 
neglected during the first part of the dissolving process.) The propor- 
tionality constant depends on the efficiency of the agitator. With a 
simple paddle wheel stirrer the following run was made: 

Volume of solution = 4000 cc. Weight of crystals K 2 Cr 2 O 7 = 1000 
gms. Temperature = 23 C. 

= 5 10 minutes 

c = 1.60 1.93 equivalents per litre. 

A second run under similar initial conditions with propeller and draft 
tube: 

= 7.5 

c = 1.95 

Find the ratio of the proportionality constants in the two cases. 
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27. A crystal of K 2 Cr2O 7 falls through a column of pure water. 
After falling 10 ft. it is dissolved. How far must it fall to completely 
disappear? Assume the diffusion coefficient proportional to the 
velocity of the crystal through the liquid, and the stream lines of liquid 
about the crystal smooth. 

28. A mass of crystals of potassium bichromate of substantially 
uniform size is suspended in water by continued, constant agitation. 
The amount of water used is many times that required to give a satu- 
rated solution. At the end of 10 miris. of agitation, titration of the 
solution shows that half of the K 2 Cr 2 O7 has dissolved. How much 
longer must the agitation be kept up to dissolve it all? 

29. One gram mol each of potassium persulphate and phosphorous 
acid were mixed, in excess of water, with 0.0569 gin. mol. hydriodic acid. 
A first reaction was 

K 2 S 2 O 3 + 2 HI = K 2 SO 4 + H 2 SO 4 + I 2 
which was accompanied by the consecutive reaction 
H 3 PO 3 + I 2 + II 2 O = H 3 PO 4 + 2 HI 

Both reactions are known to be of the second order. The constants 
were found to be ki = 0.0033 and k% = 0.079 respectively. Plot 
curves showing the amounts of iodine and of persulphate at time t 
(hours). 

30. Potassium permanganate and oxalic acid, in the proportion of 
two molecules of the former to five of the latter, were mixed with excess 
of sulphuric acid, manganese sulphate, and water. A first reaction was: 

2 KMnO 4 + 3 MnSO 4 + 2 H 2 O = 5 Mn0 2 + K 2 S0 4 + 2 H 2 S0 4 
accompanied by the consecutive reaction 

Mn0 2 + H 2 C 2 O 4 + H 2 S0 4 = MnS0 4 + 2 H 2 O + 2 C0 2 
The amount of oxalic acid was observed to change thus: 

When * = 2 3 4 5 6 7 8 9 min. 
H 2 C 2 O 4 = 51.9 42.4 35.4 29.8 25.7 22.2 19.4 17.3 

Determine the constants k\ and fa of the reactions, and then make a 
table showing the agreement of calculated with observed values for 
the oxalic acid. 

31. The corrosion of iron under water is known to be a function of 
the temperature, of the velocity of the water by the surface of the metal, 
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and proportional to the oxygen concentration in the water. Formulate 
differential equations for the rate of corrosion at constant tempera- 
ture. 

(a) in a closed vessel such that no oxygen is present other than that 
dissolved in the water at the start; 

(6) in an open vessel permitting the diffusion of oxygen into the 
water at a rate proportional to the difference between the oxygen con- 
centration in the air and that in the water; 

(c) in a pipe through which water flows at a constant rate. 

32. The equilibrium concentration of tannin on tanbark may be 
represented by the adsorption equation 



The tanbark is to be extracted with water flowing by the bark, the rate 
of extraction being proportional to the difference between the con- 
centration of the tannin in the solution in contact with the bark and 
the equilibrium concentration in the surface film of solution on the 
bark itself given by the above formula. Set up differential equations 
for the rate of extraction 

(a) in a continuous counter-current system; 

(6) in an intermittent system) where the bark is stationary and the 

liquid flows at a constant rate. 

33. Show that in a counter- current air-drier we have 

dW = kW(M - W) dx 

where W is the total water content of the stock passing a section at a 
distance x from the end where the stock enters per unit time. What 
can be said of the quantities k and M? 

34. It is known that within the temperature range from 200 C. to 
500 C. the rate at which a hydrocarbon oil, such as petroleum, cracks 
or decomposes into lighter hydrocarbons and carbon doubles approxi- 
mately for every 10 C. rise in temperature. A 1000 barrel still is 
charged with oil, brought up to 350 C., and distilled off in the presence 
of injected steam at a uniform rate over a period of sixty hours. The 
cracking is estimated at 40%. How much more rapidly would it be 
necessary to distill it at 400 C. to keep the cracking the same? 

35. With the data of Prob. 34, the oil can if desired be distilled at a 1 
progressively increasing temperature. If the distillation is started at 
300 C. and finished at 400 C.. the distillation rate being constant and 
the rate of increase of temperature being also constant, what time must 
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be taken for the distillation in order to keep the cracking the same as 
before? 

36. In commercial practice, crude petroleum is refined in a series of 
10 continuously operated stills, each holding 1000 barrels. The feed, 
10,000 barrels per day, enters one end of the first still, flows through it 
and leaves at the opposite end, whence it enters the beginning of the 
second still, and so on through the series. From each still 900 barrels 
per day is distilled off, and the final residue of 1000 barrels is discharged 
from the last still. The temperature at which the distillation starts in 
the first still is 300 C., in the second 10 C. higher, and so on. A 
uniform temperature gradient may be assumed through each still. 
Using the data of Prob. 34 so far as necessary, make a table or a graph 
showing (a) the amount of cracking realized in the stock leaving each 
successive still and (b) the time the stock leaving each still has been 
under distillation. 

37. CO is produced by passing pure C02 at atmospheric pressure 
up through a bed of coke kept by electrical heating at a substantially 
constant temperature sufficiently high so that at equilibrium the re- 
duction is practically complete. A two foot depth of fuel gives a gas 
60% CO. (a) How deep must the bed be to give 95%? (6) How 
much more slowly must the gas enter to give 95% with the same depth 
of fuel? Discuss the validity of the assumptions made in your solution. 

38. The ammonia is to be removed from a dilute aqueous solution 
at 20, containing one per cent ammonia, by aeration. If 99% of the 
ammonia is removed in this way how many cu. ft. of air must be cm- 
ployed per pound of ammonia removed? Given that the vapor pressure 
of water at 20 = 17.4 mm. Hg. and v.p. of NH 3 above a 1% solution 
at 20 = 16 mm. Hg. Assume that v.p. of NH 3 is proportional to the 
concentration of the solution. 

39. Repeat Prob. 38, supposing the original solution to be 5%. 

40. Sugar is fermented in open vats to alcohol at 25 C., the original 
sugar solution being of such strength as to produce a 10% product. 
At the temperature of the process the product exerts a pressure of 32 
mm. of which 11.83 is alcohol and 20.17 is water, the composition of 
vapor over a 10% alcohol solution being 60% alcohol by weight. One 
mol of CO 2 is formed for every mol of alcohol. Neglecting any water 
loss, estimate the percentage of the total alcohol formed which is lost 
during fermentation. 

(a) assuming that the partial pressures are not lowered by the sugar; 

(b) assuming the partial pressure of the alcohol and of the water 
depressed by the presence of the sugar in accordance with Raoult's 
law. 
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41. A dust particle settles in air at a substantially constant rate. 
The rate is a function of the size of the particle, the smallest falling at 
the lowest rate. One of the commercial processes for the removal of 
dust from air is by means of a settling chamber consisting of two hori- 
zontal parallel plates between which the air carrying dust flows hori- 
zontally at a rate slow enough to enable the particles to fall to the bot- 
tom. The velocity of the air along the walls is zero and increases toward 
the middle of the flue according to the formula from hydrodynamics 

2 _ 7/2 = 523. 

where a is half the distance between the plates, y is the distance of an 
air particle from the middle of the flue, and x is the velocity of the par- 
ticle, while b is a constant. Letting u denote the constant vertical 
velocity of the smallest dust particle, determine the length of flue needed 
to rid the air completely of dust, the constants a, 6, and u being 
known. 

42. With the tanks of Prob. 4, suppose that water flows in at the 
top of each tank at the constant rate of 4 cubic feet per minute, and 
furthermore that the second tank has an orifice in the bottom two 
inches in diameter. The first tank being initially full and the second 
empty as in Prob. 4, set up differential equations satisfied by the re- 
spective depths hi and hz at time 6. Study the changes in water level 
in the tanks. 

43. A hot fluid is carried in a horizontal direction through a rec- 
tangular chamber by a very large number of very small pipes evenly 
spaced and near together. A cold fluid passes in a vertical direction 
through the space within the chamber outside the pipes. Set up 
differential equations satisfied by the temperatures of the two fluids 
on their way through the chamber, assuming a steady state. (Processes 
in which two different gases or liquids flow at right angles are of great 
practical importance.) 

44. A coke oven consists essentially of two parallel walls of large 
extent, 18 inches apart, maintained at the constant temperature of 
2600 F. The stock when first placed in the oven undergoes rapid 
changes of a complex character, but at a time which we may call 
= it may be assumed that the temperature of the stock is substan- 
tially uniform at 300 F. Thereafter it obeys the usual laws of heat 
flow. Plot a series of curves showing the temperature gradient through 
the stock at intervals of one hour. Given that the effective conductivity 
is 0.1 in B.t.u., the unit of length being one foot and the unit of time 
one hour. The density is 100 Ibs./cu. ft. and the specific heat is 
about 0.5. 



MISCELLANEOUS PROBLEMS 105 

45. The specific heat of graphite is as follows: 

TC. -244 -186 -50.3 +10.8 +138.5 

Sp. heat 0.005 .027 0.1138 0.1604 0.2542 

TC. +249.3 +641.9 +977.0 

Sp. heat 0.3250 0.4450 0.4670 

By graphic integration plot a curve showing the total heat content 
of graphite as a function of temperature. 
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Oil wells, 99 
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of a process, 28 
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Outflow, formula for, 15 



Partial differential equations, 77 
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intermittent, 70 

of first order, 9, 13 

of second order, 27 
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Sulphur extraction, 29 
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Thorium X, 26 
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Volatile liquids, 87 Work, minimum, 90 
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